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CHAPTER  I 


INTRODUCTION 

The  Introduction  of  solid  state  arrays  will  have  a major 
impact  on  at  least  three  areas  of  semiconductor  electronics: 

1.  As  an  analogue  shift  register  or  delay  line 

2.  As  a serial  memory  for  binary  data  stiorage 

3.  As  a solid  state  imaging  device. 

Imaging  arrays  have  attractive  military  applications  to  terminal 
missile  guidance,  fire  control.  Image  reproduction,  target  or 
pattern  recognition  and  similar  areas. 

Solid  state  devices  (SSlDs)  are  attractive  in  many  areas 
because  arrays  of  charge  coupled  elements  can  be  operated  without 
the  electron  be.im  scanning  mechanism  and  high-voltage  vacuum  tech- 
nology of  conventional  vldicons.  Most  devices  are  coa^>atlble  with 
other  MOS  devices  or  bipolar  technology  which  might  be  used  as 
peripheral  circuitry,  whether  it  be  signal  processing,  clocking 
waveshapes,  external  amplifiers,  or  encoders. 

One  popular  Imaging  device  which  shows  much  promise  in  state- 
of-the-art  lamglng  is  the  charge  coupled  device  (CCD).  The  CCD 
operates  by  c mechanism  of  charge  storage  and  transfer  of  these 
packets  of  minority  charge,  which  represent  analogue  or  semi-analogue 
signals  from  under  an  array  of  HOS  control  gates.  Information  in 
this  form  is  then  transferred  along  the  silicon  surface  in  clocked 
shift  register  fashion  by  sequential  manipulition  of  the  voltages  on 
the  control  gates.  This  transfer  takes  place  in  the  bulk  silicon 


Just  beneath  the  silicon  surface.  There  are  basically  two  types  of 
CO),  surface  channel  (SCCD)  and  buried  channel  (BCCD) . In  either 
case,  the  minimun  potential  energy  of  the  depleted  wells  is  deter- 
mined by  tjie  voltage  applied  to  the  gate  electrode,  so  that  the 
appropriate  manipulation  of  the  phase  voltages  causes  the  charge 
packet  transfer  along  the  line,  since  the  charge  always  moves  to  the 
local  potential  miniisuai.  One  section  of  a three-phase  CCD  is  shown 
in  Figure  1 to  illustrate  the  potential  well  and  development  of  the 
charge  packet  under  the  control  electrode. 


Figure  1.  Section  of  a Three-phase  CCD. 
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In  the  figure.  V represents  the  electrode  gate  voltage.  Once 
8 

a charge  packet  has  been  produced  on  the  surface  of  the  device,  It 
■ay  be  transferred  along  the  structure  as  Illustrated  In  Figure  2. 

At  tlae  t^,  a charge  packet  Is  developed  under  the  electrode  con- 
trolled by  vhlch  has  an  applied  voltage  of  • At  this  saae 

tlae,  electrodes  and  <^3  are  aalntalned  at  the  resting  potential, 

ss 

At  tlae  tj,  the  electrode  controlled  by  Is  pulsed  to  , 

cc 

producing  a potential  well  under  the  control  electrode.  Since  the 
electrodes  are  claaely  spaced,  the  charge  begins  to  aove  slowly  froa 
the  ^2  electrode  to  the  controlled  electrode.  The  voltage  on  the 
^2  electrode  la  then  veduced  to  with  a slowly  falling  fuige. 

This  slowly  falling  edge  provides  the  finite  amount  of  tlae  necessary 
for  the  charge  carriers  zo  diffuse  across  the  width  of  the  electrode. 
Figure  2(c)  shows  the  charge  transfer  complete  at  tlae  t2  * with  the 
charge  stored  In  the  %>ell,  or  depletion  region,  under  the  electrode 
controlled  by  aust  be  kept  at  a relatively  low  potential  for 

this  entire  tlae  to  prevent  the  backflow  of  charge.  Thus,  In  a 
three-phase  device,  three  array  elements,  or  electrodes,  are  neces- 
sary to  store  and  transfer  one  charge  packet.  For  this  reason, 
these  three  control  electrodes  are  usually  referred  to  as  "one 
element"  of  the  CCD.  This  will  be  the  notation  used  throughout  the 
body  of  this  work. 

Thus,  the  electrodes  of  the  three-phase  CCD  array  are  arranged 
in  triplets  and  are  connected  sequentially  to  the  drive  lines  carry- 
ing the  three-phAse  clocking  waveforas,  or  driving  pulses.  These 
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complete  clocking  pulses  are  Illustrated  in  Figure  3.  It  should  be 
noted  that  there  can  be  charge  under  every  third  electrode  in  an 
array  and  that  application  of  the  clocking  waveforms  causes  these 
charges  to  move  simultaneously.* 

The  theoretical  aspects  of  charge  coupled  device  operation  will 
not  be  dealt  with  explicitly,  as  they  provide  little  insight  into  the 
operation  parameters  of  the  device  in  the  manner  of  which  they  were 
examined.  However,  it  is  helpful  to  observe  how  the  potential  well 
depth  is  related  to  the  applied  gate  voltage,  and  how  the  charge 
transfer  time  is  related  to  element  dimensions.  The  potent!^  well 
depth  may  be  expressed  as 


where 


B 


IT- 


Qss 

^ox 


ox 


(2) 

(3) 


■ gate  voltage 

Q * oxide  charge  per  unit  area 
88 

* oxide  capacitance  per  unit  area 
4 - electronic  charge  In  coulombs 

w 3 

■>  substrate  doping  level  in  cm 

e » silicon  dielectric  constant  in  F/cm 
s 

V " potential  well  depth 

O 
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C ■ oxide  dielectric  constant  in  F/cm 
ox 

X o:i'lde  thickness  in  cm. 

ox 


The  -well  depth  is  then  determined  by 


AV, 


c + c. 

ox  d 


(4) 


where  if  the  depletion-layer  canacltance,  and  is  the  signal 

charge . 

The  most  important  figure  of  merit  for  a charge  transfer  device 
is  the  trancfer  efficiency,  n , the  fractional  part  of  charge  that 
is  transferred  from  one  electrode  to  the  next.  Likewise ,"^^the  portion 
which  is  not  transferred,  the  transfer  Inefficiency,  or  transfer 
loss,  is  denoted  by  e . Thus, 


n + e - 1 


(5) 


Since  CCD  imaging  devices  require  large  numbers  of  transfers, 
e must  be  very  small  to  prevent  excessive  signal  degradation. 

The  transfer  efficiency  is  dependent  on  two  aspects  of  the 
device — the  transfer  time  and  trapping  effects.  The  trapping  effects 
con  be  minimized  by  the  continuous  circulation  of  a small  amount  of 
charge  ("lOZ  of  full  signal  level).  This  background  charge  is 
generally  called  "fat  zero",  and  helps  to  keep  the  fast  states  con- 
tinuously filled  so  that  no  states  are  empty  to  trap  charge  when  a 
full  well  signal  arrives. 

The  transfer  time  is  essentially  determined  by  either  thermal 
diffusion  or  fringing  field  drift.  Both  of  these  mechanisms  cause 
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au  exponential  decay  of  charge.  The  time  constants  for  this  decay 
determine  the  efficiency  that  can  be  achieved  at  any  particular 
clock  frequency,  since  this  clock  frequency  determines  the  time 
available  for  charge  transit. 

The  time  constant  for  thermal  diffusion  can  be  determined  by 
physical  properties  of  the  gate.  The  motivation  Is  provided  by  the 
thermal  voltage  kt/q  applied  acroas  the  gate  length.  The  transit 
time  is  therefore 


T 


th 


yE 


L 

kt  / , 
li  /qL 


(6) 


To  achieve  efficiencies  of  99.99%  (e  - lO”**)  rer(,ulre8  about  tan  time 
constants  (lO'**  - e~^’) . Thus  for  typical  values  (L  • lOum  and 
D - lOcm^/sec),  efficiencies  of  99.99%  ccn  be  achieved  at  approxi- 
mately 1 MHz.  clock  rates,  assuming  thermal  diffusion  Is  the  only 
mechanism  responsible  for  charge  transfer.  This  speed  is  Insuffic- 
ient for  many  applications. 

The  transfer  of  charge  Is  also  aided  by  the  externally  applied 
voltages  In  the  form  of  Induced  fields.  The  transit  time  constant 
for  frlnglng-fleld  drift  Is  then  given  by 


T 


f 


L* 

3.2  yVx 


OX 


(7) 


where  E Is  the  minimum  frlnglng-fleld  possible,  given  by 
n 


» V *ox 

m L L 


(8) 
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V Is  the  applied  voltage,  is  the  oxide  thickness,  and  L Is  the 
gate  length.  Again,  for  10  pm  gates,  the  thermal  time  constant  Is 

10  ^ seconds.  A8suuln^.^y  - 400,  * 1000  A , and  V ■ 10  volts,  the 

fringing  field  time  constant  In  7.8  x 10  seconds,  or  a factor  of 
13  lover.  Thus,  for  10  pm  gate  lengths  and  low  substrate  doping, 
fringing  field  drift  results  In  about  10-15  MHz.  operation.  Since 
Is  dependent  on  L^,  the  increeae  in  speed  over  that  possible  vitVi 
thermal  diffusion  Increases  as  gate  length  Is  decreased.^ 

The  maximum  amount  of  signal  charge,  Q,  which  can  be  stored 
under  any  one  control  electrode  and  transferred  within  the  CCD  Is 
given  by: 

Q - kCV 

where  C Is  the  total  oxide  capacitance  of  the  storage  electrode  and 

V Is  the  voltage  swing  on  the  control  electrode.  At  full  well,  a 
constant  fraction,  k,  of  this  is  minority  charge,  where  k»  0.5  for 
typical  surface  channel  structures.^ 

The  transfer  efficiency  can  be  calculated  from  the  relative 
magnitude  of  the  residual  "trailing"  charge  Aq  to  the  full  charge 
packet  Q : 

^ - ne  (9) 

where  n Is  the  number  of  transfers.  For  example,  for  surface 
channel  devices  e Is  typically  in  the  range  5 x 10“’  to  lO''"  (for 
medium  frequency  operation)*  The  values  for  the  ri£  product  can  be 


i 
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shown  to  be  necessarily  less  than  0.1  for  maximum  bandwidth  capa- 
hllltles  of  the  device  to  be  realized. 

For  Imaging,  an  optical  system  focuses  the  light  onto  the  CCD. 
The  incident  light  quanta  enter  the  substrate  and  impart  their 
energy  to  the  silicon  causing  the  generation  of  free  charge  carriers. 
Depending  on  device  structure,  the  quanv*:a  may  enter  between  the 
control  electrodes,  through  transparent  electrodes  or,  in  specially 
thinned  devices,  through  the  back  face  of  the  substrate.  The 
generated  carriers  collect  as  a charge  pattern  under  the  array  of 
electrodes.  This  pattern  is  a semi-analogue  replica  of  the  radiation 
of  light  intensity  across  the  original  image.  Seml-cnalogue  is  a 
term  that  will  be  used  to  describe  the  sampled  analogue  image  stored 
in  the  device.  The  charge  pattern  ciin  be  extracted  from  the  CCD 
using  the  previously  described  sequential  clocking  technique  by  one 
of  three  different  types  of  transfer.  It  appears  as  a train  of 
pulses  whose  amplitudes  vary  with  the  grey  scale  of  the  image.  The 
array  is  thus  self-scanning. ^ 

Three  options  of  device  organization  are  possible.  These  are 
frame  transfer  (FT),  interline  transfer  (ILT),  and  line  transfer  (LT). 
These  three  options  are  illustrated  in  Figure  4.  A listing  of  these 
options  end  design  icisues  are  also  given  in  Table  1. 

Frame  transfer,  which  is  a low  frequency  continuous  transfer, 
and  line  step,  which  is  also  low  frequency  but  discontinuous,  suke 
up  the  vertical  transfer  processes. 


I 

i 

I 

I 

i 

1 

I 

j 

I 

s 

I 

I 


1 

j 


The  frame  transfer  orgs"'*  sation  requires  r ubstantlally  mure 
total  device  ares  to  accomsKMlste  s buffer  store,  which  is  a major 
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Figure  4.  Three  Basic  Types  of  Organization  for  a 
Charge-Coupled  Area  Image  Sensor. 


Table  1.  Design  Issues  and  Options  for 
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disadvantage  compared  to  the  other  two  devices.  Hov^ever,  It  Is 
functionally  the  simplest.  The  line  transfer  organisation  employs 
a scan  generator.  The  Interline  transfer  structure  requires  separate 
photosensitive  elements  and  shift  registers.  The  ILT  simplicity  is  a 
definite  advantage  with  regard  to  interlace  and  vertical  resolution.** 

In  the  frame-transfer  structure,  the  top  half  of  the  chip  Is 
photosensitive.  Field  A Is  formed  by  collecting  photoelectrons  under 
the  electrodes  for  approximately  1/60  second.  This  charge  con- 
figuration la  shifted  into  the  shielded  storage  register  In  typically 
1/600  second.  Field  A Is  then  read  out  a lir^e  at  a time  while  Field 
B Is  being  formed  by  collecting  photoelectrons  under  the  ^2 
trodes. 

In  the  Interline  transfer  structure,  the  shielded  vertical 

readout  registers  are  Interlaced  with  the  photosensitive  lines. 

Because  the  Integrating  cells  and  shift-cut  cells  are  separate,  the 

effective  Integration  time  for  both  Field  A and  Field  B Is  1/30 

second.  After  collecting  photoelectrons  In  Field  A for  1/30  second, 

the  charge  configuration  Is  shifted  Into  the  shielded  registers  and 

down,  a line  at  a time,  into  the  horizontal  output  register.  When 

Field  A is  co8q)letely  read  out  <1/60  second) , Field  B Is  shifted  Into 

the  shielded  registers  and  out.  It  is  important  to  n.)te  that  the 

effective  integration  time  for  interline  transfer  structure  is  twice 

s 

that  of  the  frame-transfer  structure  or  line  transfer  structure. 

Research  has  been  perforswd  by  several  aerospace  companies  to 
develop  trackers  which  utilize  imaging  arrays  such  as  the  CCD  as 
sensors.  The  Avionics  Laboratory,  Wrlght-Patnerson  Air  Force  Base, 
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Ohio  initiated  a proijran  with  the  purpose  of  determining  the  per- 
formance characteristics  of  array  trackers,  providing  useful  test 
data  to  designers  and  exploring  the  potential  of  this  new  technology. 

Tests  are  performed  on  array  trackers  from  three  companies: 
Honeywell  Incorporated,  McDonnell  Douglas  Astronautics,  and  Martin 
Mariett.:  Aerospace. 

The  testa  are  performed  using  a Zoom  Optical  Target  Simulator 
(ZOTS) , an  Irstrument  designed  specifically  for  testing  electro- 
optical  tracking  syatems.  The  ZOTS  provides  simulated  parameters, 
such  as  brightness,  contrast,  motion,  and  range  closures.  Motion  of 
the  target  carriage  provided  in  x,  y,  and  diagonal  direction  by  a 
servo  control  system. 

Both  the  ZOTS  and  the  tracker  under  test  are  interfaced  into  a 
Data  Aquisitlom  System,  and  recorded  in  digital  form  on  taagnetic  tape. 
The  tape  format:  is  9'-track,  800  bpi,  and  IBM  compatible.  The  run 
number,  run  titae,  x and  y tracker  outputs,  ZOTS  analogue  voltages, 
atd  other  information  are  fecorded  for  each  tracker  frame. 

These  data  are  used  in  portions  of  the  parameter  investigation  for 
CCD,  charge  injection  device  (CID),  and  photodiode  arrays  utilized  in 
target  and  pattern  tracking. 
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CHAPTER  II 

PROBLEM  DESCRIPTION 


The  problea  is  three-fold.  There  is  an  analysis  of  data 
obtained  from  the  available  trackers,  a parameter  Investigation 
utilizing  these  data  as  well  as  an  extensive  literature  search,  and  a 
simulation  of  the  CCD  utilizing  parameter  variations.  As  an  exten- 
sion of  this  work  the  vxdeo  data  are  encoded  In  order  to  reduce 
redundancy  Incurred.  The  Huffman  coding  of  run-length  sequences 
accomplishes  this  task  easily. 

The  analysis  task  consists  of  Identifying  the  sources  of  tracker 
error  and  where  possible  eliminating  the  systematic  errors.  A 
thorough  Investigation  of  the  tracker  responses  to  variations  In  the 
parameters  In  question  Is  accoi^>ll8hed  as  well  as  the  removal  of 
system  errors,  revealing  the  actual  tracker  error. 

Since  the  data  are  stored  on  9-track  magnetic  tape,  a program 
written  to  retrieve  the  data  and  transform  the  information  Into  a 
usable  form  la  utilized.  The  words  ace  converted  to  a usable  length 
and  manipulated  on  the  available  Unlvac  1106  Multiprocessor  Facility. 
A major  portion  of  this  work  is,  however,  dedicated  to  the  Identi- 
fication of  these  error  sources  and  systematic  errer  removal. 

The  results  are  used  as  a partial  source  of  parameter  study 
for  the  second  portion  of  the  study. 

A synopsis  of  the  computer  software  used  in  the  error  analysis 
is  available  in  Appendix  B of  this  text. 
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Several  parameters  are  considered  and  trade-offs  examined  to 
evaluate  the  performance  and  establish  the  practicality  of  the  SSID 
as  a tracking  sensor.  There  are  several  optical  characteristics  of 

the  device  itself  that  are  significant.  For  example,  polyslllcon, 

% 

which  is  used  extensively  to  form  the  optically  transparent  electrodes, 
does  provide  some  optical  attenuation.  Also,  when  multilayers  of 
polyslllcon  electrodes  are  employed,  as  is  the  case  in  the  manu- 
facture of  some  devices,  it  is  possible  that  the  layers  may  act  as 
an  undesirable  filter. 

From  a systems  standpoint,  there  are  a number  of  considerat 'Lons 
that  affect  the  eventual  design  decision.  The  close  coupling  of  the 
photo  elements,  which  is  vital  for  the  required  charge  transfer,  to 
an  output  detector  creates  the  atmosphere  for  a blooming  problem. 

A heavy  optical  overload  results  in  bloosting;  in  fact,  a small, 
intense  spot  can  create  a disturbance  tia  extensive  as  to  cover  the 
entire  sensor.  Partial  cures  i'or  this  blooming  are  proposed  and 
evaluated  for  the  tracking  or  imaging  application. 

A consideration  along  this  same  line  is  that  ot  the  intensity 
of  the  target  and  the  background  field  of  view.  Besides  the  danger 
of  high  intensities  causing  blooming,  there  is  the  lover  limit  on 
target  Intensity  which  permits  detection  by  the  sensor  or  tracking 
algorithm.  This  means  that  there  must  be  a significant  difference 
betwen  the  iat‘«nslty  of  the  target  and  the  intensity  of  both  the 
background  and  the  nolae  Inherent  in  the  device  Itself. 

Inherent  noise  can  dagrade  system  performance  while  tracking. 

It  develops  in  the  following  manner.  If  the  object  moves  from  one 
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active  elenent  or  group  of  active  elements  to  another , and  if  the 
amount  of  charge  generated  and  stored  in  the  second  element  or  group 
is  not  comparable  to  tliat  of  the  preceding  element  or  group,  the 
sensor  might  misjudge  the  direction  of  travel  of  the  target,  or 
lose  the  target  altogether.  Nonunlfomltles  of  this  type  are 
usually  caused  by  limitations  in  the  manufacturing  process  for  the 
device. 

The  largest  single  cause  of  nonuniformity  In  sensor  output  Is 
due  to  material  properties.  The  generation  of  photo  electrons 
varies  according  to  the  homogeneity  of  the  bulk  material  In  which 
the  photons  are  absorbed.  As  devices  get  smaller  and  more  dense 
the  photolithographic  process  used  In  fabrication  becomes  a limiting 
factor.  Irregularities  become  dependent  on  the  irregularity  in  the 
geometry  of  the  sensing  cells.  Many  companies  have  thus  gone  to 
electronbeam  lithography  to  provide  smaller  dimensions  In  cell  size 
and  Intercell  spacing. 

The  amount  of  nonunlformlty  which  can  be  tolerated  is,  of 
course,  dependent  on  device  utilization.  In  tracking,  as  will  be 
shown  later,  the  nonuniformities  only  becosie  a valid  concern  when 
contrasts  and  Irradlances  approach  extremes.  If  Inherent  noise  Is 
nresent,  schemes  must  be  employed  to  reduce  Its  effect  on  the  video 
output  of  the  device.  For  instance,  corrections  can  be  Incorporated 
In  the  tracking  algorithm  or  In  the  external  circuitry. 

This  Inherent  noise  Is  compounded  by  the  Introduction  of  a 
fixed-pattern  "noise"  on  the  signal.  This  fixed-pattern  noise 
Is  imposed  by  an  effect  known  as  dark  current  nonunlformlty. 
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This  dark  current  Is  thermally  generated  in  the  silicon  and  also 
limits  the  storage  time  of  the  CCD.  'The  limitation  on  storage  time 
determines  the  lower  limit  on  scan  rate  of  the  sensor.  This  storage 
time  is  very  short  at  room  temperature,  but  cooled  devices  are 
capable  of  storage  times  of  up  to  an  hour  or  two. 

Tl<e  scan  rate  of  the  sensor  determines  the  maximum  imape  velocity 
of  a trackable  target.  The  scan  rate  (data  update)  must  be  somewhat 
faster  than  the  estimated  velocity  of  the  target  image  across  the 
face  of  the  device.  If  not,  the  target  will  be  lost  between  frames. 
On  the  upper  end  of  the  scan  rate,  limitations  are  imposed  by  the 
integration  time  of  the  SSID.  Integration  time  is  a fundamental 
property  of  the  device.  If  the  device  is  scanned  faster  than  the 
integration  tine  necessary  to  adequately  describe  the  scene  In  quan- 
tities of  charge,  blurring  will  result,  and  target  definition  suffers. 

As  scan  rates  and  frame  rates  are  increased,  two  interrelated 
quantities  becoise  crucial  to  system  performance.  These  are  transfer 
time  and  transfer  efficiency.  The  loss  of  the  small  portione  of 
information  and  the  subsequent  addition  of  these  residual  charges  to 
trailing  elements  tend  to  degrade  video  output.  Therefore,  scan 
rates  must  allow  for  adequate  charge  transfer.  System-wise,  higher 
scan  rates  need  wider  bandwidth  and  wide  band  noise  may  prove 
undesirable. 

Another  area  of  concern  is  target  acquisition.  This  aspect  is  a 
characteristic  of  both  the  tracking  algorithm  and  the  device  parami^ 
ters.  Thus,  various  aspects  are  Interdependent.  Tlie  sensitivity 


of  the  SSID  to  changes  In  the  background,  or  the  amount  of  grey 
level  separation  which  is  attainable,  determine  how  well  the  target 
can  be  defined  and  tracked. 

Several  of  the  above  listed  problems  are  interrelated,  and  thus 
trade-offs  are  considered.  For  instance,  integration  time,  scan 
rate,  target  velocity,  transfer  time,  and  transfer  efficiency  are  ell 
interdependent.  Likewise,  so  are  Intensity,  blur,  sensitivity. 
Inherent  noise,  and  blooming  effects. 

It  is  necessary  to  study  the  extent  of  each  probliam  as  it  is 
introduced  by  the  device  parameters  and  as  it  affects  system  per- 
formance. These  results  can  be  used  to  determine  corrective  measures. 
The  data  described  previously  are  used  as  an  experimental  basis  or 
applied  starting  point  for  the  system  effects  of  device  parameters. 
They  are  used  for  analysis  and  Investigation  of  projected  system 
performance. 

Portions  of  the  projected  performance  are  attained  from  a com- 
puter simulation  of  the  CCD  tracker  array.  The  fortran  program 
utlllz'^8  transfer  efficiencies,  target  contrasts,  background  elimi- 
nation, target  pattern  acquisition  and  grey  level  adjustment  in  each 
array  element.  The  simulation  provides  Information  as  to  the 
degradation  of  the  picture  Image  as  the  analogue  replica  of  each 
scene  la  clocked  from  the  array,  as  well  as  the  feasibility  of 
target  acquisition  on  large  arrays.  The  target  acquisitions  can 
be  utilized  as  a contrast  tracker  or  a pattern  recognition  schesie. 

The  program  Is  so  \n:ltten  that  the  Information  clocked  from 
the  array  can  be  encoded  directly  using  Huffman  encoding  to  reduce 
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scene  redundancy  aa  well  as  the  space  necessary  fo^  Inforaiatlon 
storage.  This  portion  of  the  work  Is  backed  by  previous  work  done 
In  the  area  of  redundancy  reduction  by  the  coding  of  run- length 
sequences.  Any  code  can  be  used  since  the  scheme  uses  the  table 
look-up  procedure  for  coding. 
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DATA  ANALYSIS 

Each  tracker  tested  In  this  study  is  listed  in  the  order  In  which 
the  original  tests  are  made.  The  tracking  algorithms  are  proprie- 
tary to  the  indlvldtial  con^aniest  and  hence,  are  only  discussed  in 
a limited  fashion.  The  programs  used  for  extraction  and  reduction 
of  the  data  for  this  study  are  located  in  Appendix  B.  A synopsis 
of  each  of  these  programs  is  also  given. 

TRACKER  NUMBER Honeywell  Solid  State  Imaging 
Seeker  (Photodiode) 

The  data  from  two  Honeywell  array  trackers  are  available  for 
investigation.  Both  utilize  the  same  tracking  algorithm,  but  have 
different  detector  arrays.  Tracker  Number  1 uses  a Reticon  SO  x 50 
array  of  photodiodes  in  matrix  form.  Each  cell  stores  charge  in  the 
p-n  Junction  and  is  read  serially  once  during  each  frame. 

This  device  is  designed  as  a terminal  guidance  seeker  for  a 
shoulder-launched  anti-tank  weapon.  Its  tracking  algorithm  is  a 
correlation  technique  utilizing  information  from  a specific  number 
of  array  cells  in  the  form  of  several  levels  of  "grey"  or  internal 
target  contrast. 

TRACKER  NUMBER  2.  McDonnell  Douglas  Mosaic  Tracker 

Tracker  Nuaber  2 utilizes  a General  Electric  100  x 100  charge 
injected  device  (CID)  array.  This  sensor  consists  of  a matrix  of 
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coupled  MOS  capacitors  as  does  the  CCD.  Each  cell  is  read  once  per 
frame  by  Injecting  the  charge  Into  the  silicon  substrate. 

The  tracking  algorithm  of  this  tracker  consists  of  calculating 
the  mathematical  centroid  of  the  target  from  the  classical  defi- 
nition of  the  center  of  mass  of  a uniformly  dense  target.  This 
target  tracker  uses  two-level  logic.  A present  threshold  determines 
the  presence  of  the  target  over  a particular  array  element.  Thus, 
any  cell  outputting  a 1 (as  determined  by  the  threshold)  Is  Included 
In  the  tactical  centroid  calculation. 

TRACKER  NUMBER  Martin  Marietta  Area 
Correlation  Tracker 

This  Martin  Marietta  tracker  uses  a Fairchild  201  CCD  as  a 
sensor.  It  Is  a 100  x 100  matrix  of  cells  aa  Is  used  in  the  simu- 
lation later  on  In  the  text  of  this  study. 

The  tracking  algorithm  of  this  tracker  utilizes  target  signa- 
tures characteristic  of  the  target  and  background  features  within  a 
fixed  segment  of  the  sensor  field  of  view.  The  tracker  stores  a 
reference  target  area  signature  In  memory  and  derives  tracking 
Information  from  the  comparison  of  subsequent  target  signature  with 
this  reference.  The  target  reference  can  be  updated  to  enable 
tracking  of  relatively  fixed  targets  during  range  closure.  It  Is 
this  Idea  that  la  utilized  in  the  simulation. 

TRACKER  NUMBER  A.  Martin  Marietta  Point  Tracker 

The  point  tracker  once  again  uses  the  CCD  array  as  does  tracker 
number  3«  the  difference  in  the  two  being  In  the  tracker  electrorlca. 
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The  point  tracker  algorltho  processes  sensor  Infonoatlon  to  dis- 
tinguish a target  from  Its  background  and  then  calculates  the 
relative  location  of  the  centroid.  This  tracker  determines  the 
contrast  edges  of  the  target  and  adjusts  tracking  to  match  the 
target  window  with  the  size  of  the  target. 

TRACKER  NUMBER  5^  Honeywell  Solid  State 
Imaging  Seeker  (CID) 

This  Honeywell  device  also  uses  the  General  Electric  CID 
sensor.  However,  the  sensor  utilizes  a 1A4  x 288  cell  array.  This 
tracker  uses  the  same  tracking  algorithm  as  Tracker  Number  1,  the 
only  modifications  being  to  accotnnodate  the  array's  larger  size. 

The  dimensions  of  the  target  Image  on  the  sensor,  the  angular 
size  and  the  approximate  nuinber  of  sensor  cells  subtended  by  the 
target  are  presented  in  Table  2. 


Table  2 

Geometric  Target  Sizes 


Tracker  Number 

Target  Size 

Mils  Mllllradlans 

Cells  Subtended 

1 

12  X 12 

2x2 

3x3 

2 

16. A X 16.4 

2.7  X 2.7 

5 X 6.5 

3 

24.1  X 24.1 

4x4 

15  X 20 

A 

13.3  X 13.3 

2.2  X 2.2 

8 X 11 

5 

7.2  X 4.2 

1.7  X 1 

3x3 
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In  order  to  present  the  parasieter  study  perfortoed  on  the  three 
companies'  trackers • an  Identification  and  analysis  of  tracker 
error  Is  necessary. 

The  error  analysis  Is  begun  by  reviewing  saitiple  runs  and  obtain- 
ing plots  of  the  ZOTS  position  versus  time  for  Tracker  Nuiii>er  1.  The 
ZOTS  position  Is  the  true  position  of  the  target  and  Is  measured  In 
mils.  The  time  is  given  in  seconds.  The  k and  y position  on  the 
array  sensor  are  recorded  as  the  x and  y analogue  outputs.  The 
computer  program  TAPEREAD.DRIVEI  calculates  the  sensor  error  dif- 
ference E as 
8 

E - Z - T (inches)  » (11) 

8 

where  Z is  the  ZOTS  Irdlcation  of  "true"  target  position  and  T is 
the  tracker  realization  of  target  position. 

The  calculated  error  for  the  Honeywell  tracker  is  plotted  for 
investigation  in  Figure  5.  Several  tynes  of  error  are  evident  from 
the  plot.  The  error  function  has  an  average  value  which  takes  on 
a positive  value  during  the  first  pass  of  the  target  on  the  image 
surface.  When  the  target  changes  directions  and  begins  its  reverse 
passy  the  value  drops  to  a alightly  negative  quantity.  This  error, 
it  seems,  is  due  to  tracker  lag.  This  tracker  lag  might  not  be  a 
characteristic  of  the  imaging  device  Itself,  but  an  inconvenience 
encountered  in  the  use  of  certain  tracking  algorithms.  The  lag 
error  is  evident  in  all  of  the  trackers,  but  varies  in  its  extent. 

The  higher  frequency  oscillations  represent  "electronic  noise". 
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proportional  to  the  location  of  the  target.  Superimposed  on  theae 
d.  c.  voltages  is  a random  “electronic  noise*',  which  is  Interpreted 
as  “true"  target  motion  by  the  tracker.  This  nolae  is  easily  re- 
moved in  the  computer  analysis  by  simply  obtaining  a linear  least- 
squares  fit  to  the  ZOTS  output  voltage.  The  removal  of  this  noise 
reveals  the  actual  tracker  error  as  pictured  In  Figure  6.  These 
results  correspond  much  more  closely  to  the  reaction  of  an  “Ideal" 
tracker.  An  “ideal"  tracker  error  occurs  only  between  target  posi- 
tion decisions.  The  difference  error  Increases  from  zero  until  the 
tracking  algorithm  makes  the  decision  that  the  target  has  Indeed 
moved  to  another  cell  space  or  active  element.  At  the  predetermined 
time,  the  tracker  output  indicates  a move  and  the  error  drops  to 
near  zero,  where  the  procedure  is  Itterated;  resulting  in  a sawtooth 
waveform  whose  frequency  is  given  by 

similar  to  that  of  Figure  7.  v^  is  the  target  velocity  and  d is 
the  width  of  an  active  element  of  the  sensor.  There  is  strong  indi- 
cation that  this  is  Indeed  the  case,  in  part,  showing  most  of  the 
error  to  be  algorithm  oriented  rather  than  sensor  oriented. 

The  other  error  element  encountered  in  the  analysis  is  that  of 
calibration  factor  selection*  The  calibration  factors  for  scaling 
outputs  from  the  ZOTS  and  tracker  determines  the  amplitude  of  the 
sweep  for  each.  Lack  of  precision  produces  a large  triangular 
component  in  the  error  difference  as  shown  in  Figure  8.  This  same 
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Figure  7.  Explanation  of  Actual  Error  (Sawtooth). 


Figure  8.  Explanatioa  of  Calibration  Error  Versus  Time 


type  of  error  would  be  evident  if  a lag  In  response  of  the  target 
to  ZOTS  positioner,  since  the  target  would  never  cross  the  full 
expanse  of  the  face  of  the  laager.  Each  case  Is,  however,  a 
limitation  of  the  test  equlpoent  and  not  of  the  tracker  Itself.  I 
Is  not  considered  a restriction  on  the  tracker. 

Sample  runs  for  obtaining  the  tracker  error  are  shown  In 
Figures  9 through  14.  These  rune  provide  Illustrations  of  each 

type  of  error  encountered  In  the  tracker  parameter  analysis.  The 

{ ■* 

I Fortran  program  also  provides  numerical  values  for  the  error  dlf- 

i 

[ ferances,  but  an  extensive  listing  would  be  cumbersome  and  provide 

I 

I little  extra  Insight  Into  the  abilities  of  the  tracker. 


Tine  (aisec) 

>r  Honeywell  Tracker;  Irradlance  > 21.0,  Contrast 
Velocity  - 13.22,  F.  R.  « 92.6;  Run  lU. 


T 


Honeywell  Tracker;  Irr.  * 21.0,  Contr 


Tracker  Error  for  Martin  Marietta  Correlation  Tracker. 


CHAPTER  IV 


RESULTS  OF  PaRAM^ER  STUDY 

The  sources  of  the  tracker  test  and  analysis  are  a statistical 
data  reduction  program  written  for  the  available  Unlvac  1106  Multi- 
processor facility  and  ZOTS  data  from  the  laboratory  evaluation  of 
the  five  solid  state  Imaging  arrays.  Some  of  the  latter  Is  available 
from  previous  reports  and  repeated  on  the  local  facility  for  com- 
parison purposes.  The  minimum  and  maximum  values  of  Eg  are  used  to 
determine  If  loss  of  track  occurs  In  each  run.  Loss  of  track  Is 
defined  to  have  occurred  when  an  extreme  value  of  Eg  Is  significantly 
greater  than  one-half  the  appropriate  x or  y dimension  of  the  target 
Image  on  the  sensor.  Since  the  trackers  have  no  provisions  for 
rcacqulsltlon  of  the  target,  such  occurrences  are  also  defined  ns 
loss  of  track. 

The  data  shown  In  Tables  3 through  7 provide  an  analysis  of 
tracker  accuracy  of  the  five  trackers  as  a function  of  Irradlance  of 
the  tracker  lens.  These  data  are  used  to  determine  the  Irradlance 
range  over  which  each  tracker  can  successfully  track.  The  upper 
Irradlance  level  is  chosen  safely  below  the  level  necessary  to 
saturate  the  sensor  or  cause  blooming.  The  lowest  reported  Irradlance 
level  is  defined  as  the  lower  limit,  which  is  Immediately  above  the 
non-tracking  level,  x and  y motion  refer  to  Independent  target 
motion  along  that  respective  axis.  The  quantity  A Is  the  difference 
Eg  and  O Is  the  standard  deviation,  both  of  which  are  averaged  over 
all  runs  of  a particular  test  series,  expressed  In  mlllltadlans . 


Table  4.  Tracking  Accuracy  Versus  Irradlance,  Tracker  No. 
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Tracking  accuracies  as  a function  of  contrast  are  shown  In 
Tables  8 through  12.  The  contrasts  used  are  defined  In  the  fol- 
lowing manner. 


C - 100 


(13) 


where  C is  the  contrast  In  percent,  L is  the  measured  luminance  in 
foot-lamberts,  and  subscripts  T and  B refer  to  target  and  background 
respectively. 

This  test  determines  the  lowest  possible  photometric  contrast 
at  which  each  tracker  could  successfully  track. 

Tables  13  through  20  present  tracking  accuracies  as  a function 
of  target  rate  with  both  Irradlance  and  target  contrast  variable. 

The  purpose  of  this  test  Is  to  determine  the  tsaxlmum  tracking  rate 
capability  of  each  of  the  trackers  at  the  extreme  limits  of  the 
Irradlance  and  contrast  tests.  The  Initial  target  velocity  In  array 
cells/second  Is  approximately  60Z  of  the  tracker  frame  rate.  At  each 
Irradlance  or  contrast  level,  the  target  rate  Is  Increased  until 
either  the  tracker  ceases  to  track  or  the  ZOTS  limit  Is  reached. 

Sate  tests  are  performed  In  x,  y,  and  diagonal  axes  when  unrestricted 
by  ZOTS  nonuniformity  or  test  time  limitations. 

The  five  trackers  prove  to  be  very  accurate,  reliable  and  to 
have  wide  dynamic  ranges  In  both  Irradlance  and  target  contrast. 

These  capabilities  are  presented  In  Table  21.  The  Irradlance  capa- 
bilities are  expressed  as  the  dynaadc  ranges  of  successful  operation: 
since  absolute  limits  may  be  modified  by  changes  In  tracker  optica 
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Table  9.  Tradtlng  Accuracy  Versus  Ccatrast,  Tracker  Ko. 
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Table  17.  Tracking  Accuracy  Veratis  Diagonal  Target  Rate,  Tracker  No. 


Table  18.  Tracking  Accuracy  Versus  k-AxIs  Target  Rate,  Tracker  No. 
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Table  19.  Tracking  Accuracy  Versus  y-axls  Target. 
Bate,  Tracker  No.  A 

Error  In  Mllllradlans 


Ir rad lance 
yWatt/cm^ 

Contrast 

X 

16.1  nrad/sec 

A a 

28.9  mrad/sec 

A a 

0.211 

77.4 

0.046 

0.158 

0.032 

0.211 

0.0278 

82.1 

0.071 

0.167 

0.187 

0.191 

0.0075 

78.5 

0.155 

0.171 

0.134 

0.207 

0.0055* 

80.2 

0.120 

0.236 

0.103 

0.324 

0.211 

3.5 

0.101 

0.177 

0.072 

0.203 

0.214 

2.2 

-0.089 

0.208 

-0.063 

0.201 

0.218* 

0.04 

0.043 

0.164 

0.066 

0.214 

Notes:  Tracker  frane  rate,  124  fraiaes/sec 
a Adaptive  gate  disabled 
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Table  20.  Tracking  Accuracy  Versus  Target  Rate 
Tracker  No.  5 

Error  in  Milllradlans 


Irradlance 

Contrast 

Velocity 

yWatt/cm^ 

% 

Arad/ sec 

.4 

14.9  0 

.4 

15.7  0 

.4 

16.1  0 

.5 

3.48  0 

.5 

6.97  0 

.6 

4.65*^  0 

Notes: 


0.6A7 

0.537 


Tracker  frame  rate,  30  franea/sec 
Target  motion  In  x-axls 

a Lost  track  40%  of  runs,  would  not  track  at  16.6  mrad/sec 
b Lost  trc:ck  20%  of  runs,  would  not  track  at  greater 
target  velocities 

c Would  not  track  at  5.23  mrad/sec 


Table  21.  Array  Tracker  (capabilities 


Tracker  Frame  Rate 
No.  (f/sec) 


Irradlance 

«*>>»• 


Min  Contrast 
% 


Max  Target 
Rate 

mrad/sec 


1 

92.6 

10 

1.5 

- 

2 

46 

12 

12.2 

>54 

3 

119 

19 

2.1 

>36 

4 

124 

26 

3.6 

>39 

5 

30 

11 

>1.7 

16.1 
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and/or  frame  rate.  The  maxlaun  target  rate  for  trackers  2,  3,  and 
A are  not  determined  due  to  rate  limitations  of  the  ZOTS. 

One  major  advantage  of  Tracker  Number  1 is  that  it  requires  no 
adjustments  throughout  the  complete  series  of  tests  at  a given 
frame  rate.  Trackers  Number  3 and  4 require  adjustments  tor  the 
conversion  from  a high-irradiance  range  to  a low-irradlance  range. 
Trackers  Number  2 and  5 require  adjustments  for  each  parametric 
change  of  Irradiance  or  contrast. 

The  tracking  accuracy  of  a particular  tracker  (except  Tracker 
Number  2)  is  essentially  Independent  of  irradiance  until  the  marginal 
level  is  approached.  The  tracking  accuracy  is  usually  degraded  by 
10  - 20Z  at  the  margiual  level.  Tracker  Number  2 is  the  only  tracker 
exhibiting  greater  errors  than  0.36  mllllradlans  or  greater  standard 
deviation  than  0.63  milllradians  under  all  irradiance  conditions. 

The  accuracies  of  Tracker  Number  2 are  much  lower  and  decrease  with 
irradiance  level. 

The  tracking  accuracy  of  the  trackers  (excepting  Tracker 
NuBiber  2)  are  less  than  0.51  mllliradians  with  standard  deviations 
of  less  than  0.54  mllliradians  for  all  contrasts.  The  accuracy  is 
again  essentially  independent  of  target  contrast  for  contrasts  above 
the  level  at  which  the  tracker  ceases  to  operate.  Trucker  Number  2 
is  much  lower  and  decreases  with  target  contrast. 

Target  rate  effects  were  different  for  each  of  the  trackers. 


I 

1 

I 

i 

I 

i 

I 


( 


The  only  trackers  shoving  independence  to  target  ratea  are  Trackarr 
Number  2 and  5,  Trackers  Number  1 and  3 showing  linear  increases  in 
error  with  target  rate.  Tracker  Number  4 shows  a linear  increase 
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la  error  In  the  y direction • buL  tue  x-axls  error  approximately 
doubles  between  the  rates  of  8.6  and  16.1  mrad/sec.  The  x-dlrection 
error  remains  essentially  constant  for  higher  velocities.  The 
ofif^ miiiw  target  rate  capabilities  of  Trackers  Number  1 and  5 are  0.93 
times  the  theoretical  value  of  one  cell/frame,  except  for  the  tests 
at  the  marginal  irradlance  and  contrast  limits*  where  it  drops  to 
0.40  and  0.27  times  the  design  rate  limit,  respectively. 

Data  from  Tracker  Number  1 indicate  the  tracking  accuracy  at  a 
frame  rate  of  92.6  is  essentially  five  times  that  ct  23.1.  Trackers 
Number  2 and  3 demonstrate  no  significant  differences  at  different 
frame  rates.  Lower  target  velocities  account  for  the  increase  in 
accuracy  at  the  lower  frame  rate  of  Tracker  Number  1. 

At  higher  target  velocities  the  tracker  tracks  the  target  but 
lags  behind.  When  the  lag  accumulates  to  the  degree  where  the 
target  leads  by  more  than  one  call,  the  target  is  lost. 

Array  response  has  some  degree  of  noounlformlty.  This  nonuni- 
formity  contributes  to  tracker  failure  at  marginal  conditions  of 
irradlance  and  contrast.  The  array  used  in  Tracker  Number  1 has  a 
column  of  '*hot"  cells  and  the  array  of  both  Trackers  Number  3 and 
4 have  several  irregularly  shaped  areas  of  non-uniform  cells.  Under 
extremely  low  Irradlance  or  contrast  conditions,  the  target  is  lost 
if  its  image  passes  over  these  areas.  An  operational  tracker  is 
limited  by  such  nonunlformltles. 

The  five  trackers  using  the  solid  state  arrays  are  extremely 
reliable.  No  failure  of  any  of  the  tracker's  electronics  or 
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detectors  under  severel  hundred  hours  of  testing  exhibits  promising 
performance. 

Two  of  the  trackers  previously  discussed  use  a 100  x 100  CCD 
array  (Martin  Marietta).  With  the  recent  developments  in  MOS 
technology!  it  seems  apparent  that  the  charge  coupled  concept  will 
gradually  play  a more  dominant  role  in  the  imagery  technology. 

The  parameters  which  seem  to  be  sensor  oriented  rather  than 
algorithm  oriented  are  contrasts  and  transfer  time!  or  transfer 
efficiency.  The  latter  affects  the  frame  or  scan  rate  and  thus  the 
wf v4iii.^w  target  velocity.  The  dependence  of  a tracker  on  these 
parameters  depends  on  what  type  of  tracker  is  employed  and  what 
type  of  target  recognition  scheme  is  used. 

Several  systems  are  available  utilizing  acquisition  and  tracking 
algorithms!  providing  reliable  tracking  of  missile  plumes.  Here  the 
missile  plums  is  assumed  to  have  a positive  contrast  ratio  with 
respect  to  the  background.  However!  many  cases  do  not  offer  such 
conveniences.  It  is  desirable  to  track  the  target  pattern  itself » 
such  as  missile  body  trackings  aircraft  tracking!  or  projectile 
tracking.  If  a system  is  available  capable  of  such  a task!  then 
the  possibility  of  real-time  data  reduction  caa  be  considered.  The 
parameters  influencing  sensor  behavior  or  response  evidently  con- 
trol the  upper  limit  on  tracker  efficiency. 

A very  good  idea  of  how  these  parameters  Influence  device 
behavior!  is  obtained  through  Investigation  of  tha  modulation 
transfer  function.  This  modulation  transfer  function  describes  the 


59 


response  of  the  CCD  to  esch  psrsaeter  as  a function  of  frequency. 

It  is  analogous  to  the  frequency  response  of  a networic. 

Also,  aliasing  occurs  due  to  the  discrete  nature  ox  the 
sensor  eleaeuts.  In  other  words,  the  Imaging  surface  Is  not  con- 
tlnguous  but  saaq>les  the  Image  plane.  Thus,  the  spatial  frequency 
respon*.'  beyond  the  sampling  limit,  or  Nyqulst  limit,  is  a design 
Consideration.  This  characteristic,  however,  is  not  as  severe  as 
has  been  supposed.  If  It  proves  to  be  a problem,  it  can  be  removed 
by  prefiltering  or  defocusslng  the  image  to  render  the  edges  and 
lines  continuous.  This  also  tends  to  blur  the  Image  and  reduce 
edge  contrast  by  spreading  the  energy  of  the  scene  elements  to 
adjacent  areas. 

If  the  sensor  Is  to  be  used  to  detect  subresolution  sized 
objects,  such  as  point  Images  of  stellar  objects.  It  Is  a deflnlts 
disadvantage  to  have  unresponsive  Intensities  In  the  image  plana. 

It  la  preferable  then  to  have  the  sensors  contiguous,  or  for  the 
scene  to  be  prefiltered  to  that  degree.  However,  If  sensitivity  Is 
the  goal,  preflltcring  Is  not  advised.* 

The  modulation  transfer  function  (HTF)  Is  the  ratio  of  device 
output  to  Input  in  a particular  observation  domain.  In  the  frequency 
domain,  the  MIF  relates  the  amplitude  and  phase  of  a sinusoidal 
output  to  the  corresponding  sinusoidal  Input  signal.  The  argument 
gives  Che  phase  shift  and  the  modulus  gives  the  attenuation. 

Generally  the  normalization  used  Is  unity  amplitude  and  zero  argu- 
ment at  zero  frequency.  In  the  case  of  a linear  system,  the  MIF  la 
simply  the  Fourier  transform  of  the  impulse  response  function*  In 


the  case  of  a aanpled  data  ayates-f  such  a:^  a CQ>,  a single  input 
frequency  gives  rise  to  several  output  frequencies.  These  additional 
signals  are  multiples  and  sidebands  of  all  of  the  multiples  of  the 
saa^Ung  frequency,  because  the  additional  output  signals  represent 
no  additional  information  the  MTF  Is  calculated  using  the  component 
of  the  output  having  the  same  frequency  an  the  Input.  Hence  the  MTF 
becomes  simply  the  frequency  response  of  the  system.  Here,  the 
Myqulst  frequency  becomes  the  spatial  sampling  limit. 

The  MTF  Is  an  excellent  way  of  expressing  the  operational 
imaging  characteristics  of  the  CCD.  The  overall  MTF  of  the  chip  is 
composed  of  three  factors: 

1.  The  loss  of  spatial  frequency  due  to  the  geometry  of  the 
integrating  cell 

2.  The  loss  of  frequency  response  due  to  transfer  inefficiency 


3.  The  loss  of  frequency  response  due  to  the  diffusion  of 
charge  between  photon  absorption  and  photoelectron 
collection 


1- 

The  Integration  MIF  is  given  by  the  Fourier  transform  of  the 
basic  Integratloa  cell.  For  a rectangular  cell  of  length  Ax 
repeated  with  periodicity  p,  l.c.,  p is  cell  to  cell  spacing; 


the  MTF  la 
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where  is  the  Nyquist  frequency.  For  Ax  « p,  the  first  zero 
occurs  at  f ■ fj^  “ 1/p  . For  Ax  ■ p/2,  the  first  zero  occurs  at 
^vt  “ ^/P  • Figure  15  gives  the  integration  MTF  versus  normalized 
spatial  frequency.  For  the  backside  illuminated  frame  transfer 
(BIFT)  device.  Ax  * Ay  -•  p,  and  for  frontside  illuminated  interline 
transfer  (FIIT)  devices.  Ax  - Ay  - p/2, 


2.  MTF,  , 

transfer 

During  the  transfer  of  d sampled  sinusoidal  signal  along  a CCD 
shift  register,  a fraction  of  the  charge  e is  lost  from  each  of  the 
samples  at  each  transfer;  and  this  charge  is  added  to  trailing 
samples.  The  effect  of  this  dispersion  on  MTF  is  given  by 


MTF 


transfer 


exp 


j -ne[l  - CO.  j (15) 


where  n is  the  number  of  transfers. 

Figure  16  shows  the  transfer  MTF  versus  the  nv’^rmalized  spatial 
frequency  with  the  nc  product  as  the  parameter. 

The  number  of  transfers  in  the  x direction  is  the  same  for  both 
frame- transfer  and  interline-transfer  arrays.  Therefore,  v^.he  hori- 
zontal MTF  degradation  due  to  transfer  is  the  same  for  both  arrays. 
The  number  of  transfers  In  the  y direction  is  greater  for  the  frame- 
transfer  chip  by  the  amount  FH^,  where  F is  the  number  of  phases  and 
N is  the  number  of  y-dlrectlon  cells.  Therefore,  the  vertical  MTF 

y 

for  the  fraaie- transfer  chip  is  Ihwer  than  that  for  the  interline- 


transfer  chip. 


ne-0.01 
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1£  photons  are  absorbed  within  the  depletion  regions,  then  the 
collection  Is  asauoed  1002  efficient.  However,  If  photons  are 
absorbed  away  from  the  depletion  regions,  then  the  charge  configu- 
ration will  spread  as  It  diffuses  toward  the  depletion  regions  with 
a resulting  decrease  In  MTF. 

If  photons  are  absorbed  a distance  d from  the  depletion  region 
and  If  the  diffusion  length  In  the  silicon  Is  L^,  then  the  MTF  due 
to  the  diffusion  of  charge  Is  given  by 


co8h(d/Li) 

“”di£t  ■ 'HSiHTW  • 

whare  + (2llf)  . 

Figure  17  shows  the  diffusion  MTF  versus  normalized  spatial 
frequency  with  d as  the  parameter. 

The  responslvlty  of  each  photo  element  is  a description  of  the 
efficiency  of  each  cell  In  converting  Incident  photons  to  the  photo 
electrons  eventually  collected  In  the  potential  well  beneath  the 
sensor  gate. 

The  photoelemant  responslvlty  la  determined  by  the  efflcloncy 
with  which  photons  are  absorbed  and  the  resulting  photoelectrons  are 
collected.  Basically  four  mechanisms  act  to  reduce  the  photoeleaent 
r..pon.lvlty  . 

1.  Baflectlon  at  layer  irterfaces  before  the  photons  reach 


the  silicon 


0.5 


1.0 


F/F„ 

Figure  17.  Vemus  KoratellscMi  Spatial  Frequancy. 
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2.  Aborptlon  In  these  layers  before  the  photons  reach  the 
silicon. 

3.  Recombination  at  the  S1~S102  Interface  after  hole-electron 
generation. 

4.  Absorption  too  far  away  from  potential  veils  for  the 
photoelectrons  to  be  collected. 

Mechanisms  2 and  3 cause  a reduction  In  R . ..In  the  blue, 

element  ' 

mechanism  4 causes  a reduction  In  R , ..In  the  Infrared,  and 

element  ' 

mechanism  1 causes  Interference  fringes  throughout  the  spectrum. 
Mechanism  1 la  mainly  responsible  for  being  different  for 

BIFT  and  FIIT  structures. ‘ 

As  previously  mentioned,  dark  current  Is  not  completely  uni- 
form throughout  a device.  The  dark  current  signature  of  a device  Is 
obtained  by  Integrating  In  the  dark.  If  the  element- to-element  non- 
uniformity  is  , then  the  minimum  signal  that  can  be  detected 

will  be  limited  to  • ANj  . If  the  element- to-element  variation 

s d 

of  is  10  percent,  then  at  room  temperature  this  will  limit  the 
mfniimm  detectable  signal  to  about  1000  electrons.  The  dark  current 
is  a strong  function  of  temperature,  decreasing  by  a factor  of  2 
every  10*  C decrease  in  temperature. 

Three-phase  diarge  coupled  devices  with  transfer  efficiencies 
of  99.99ili  are  available  providing  operating  .frequencies  of  up  to 
10  HUs.  Since  nonunlformities  ere  diaracterlstlcs  of  the  oanufactu- 
rlng  process,  the  parameter  of  the  device  itself  in  need  of  experi- 
mental investigation  is  this  transfer  efficiency  es  it  is  related  to 
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image  degradation,  target  acquisition,  tracking  and  the  overall  MTF. 

In  the  discussion  of  the  effects  of  transfer  inefficiency  on  images, 
a single  element  with  stored  charge  Q is  observed.  This  charge  is 
placed  in  a well  (1-0),  then  after  N transfers,  the  distribution 
of  charge  in  cells  1 0,1,2,...  is  given  by 

■ r ■ 7ii-i)Tn  » - 

where  a is  the  fractional  loss  Incurred  by  the  charge  packet  in 
moving  from  one  cell  to  the  next.  The  phenomenon  is  Illustrated  in 
Figure  16  for  the  first  two  transfers.  The  fractional  loss  from  the 
charge  packet  after  N cell  transfers  ia  given  by  1 - D(N,N); 

loss  - 1 - (1  - a)“  S Na  for  Ka  « 1 (18)  j 

S 

! 

The  equation  D(M,N)  ••  D(N  - 1,N)  gives  the  value  of  N for  which  the 

first  trailing  packet  and  the  leading  packet  have  equal  amounts  of  * 

charge.  This  value  is  calculated  in  the  following  manner. 

i 

1 


D(N,N) 

- W-WIHI  « - " 

(19) 

D(N,N) 

- Hd-a)"..- 

(20) 

D(H,N) 

- (1  - a)** 

(21) 

Likewise, 


' .,*v; 

. ' + 
■-  -<■ 

• - >. 


69 


D(N-1,N) 

■ (1  - o’*-! 

(22) 

D(N-1,N) 

(23) 

D(N-1,N) 

- N(1  - a)**"^  a ; 

(24) 

thus, 

D(N,N)  - D(N-1,N)  ->  (1-a)**  - Nd-a)**"^  a . (25) 


Therefore, 


(1  - 

(26) 

! 

i 

No  - 1 - a 

(27) 

f 

1 

j 

1 

n . 1 - ct 

“ o 

(28) 

j For 

a « 1 , N - 1/a  . 

(29) 

Thus,  for  every  1/a  transfers,  the  peek  of  the  charge  distri- 
bution shifts  beck  one  cell,  1.  e. , after  2/a  transfers,  the  second 
trailing  packet  contains  the  largest  charge  quantity.  This  con- 
stitutes a delay  In  addition  to  the  delay  ties  N/f^  required  to 

c 

clock  a packet  out  of  the  CCD.** 

This  developsMnt  Is  used  to  create  a slaulatlon  of  the  CCO  as 


t 
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bit  more  complicated,  however,  since  each  cell  has  some  charge 
stored  In  its  depletion  well.  The  concept  Is  the  same,  In  order 
to  Illustrate  this  concept,  an  Image  Is  placed  on  the  100  x 100 
simulated  array.  The  test  pattern  resembles  that  shown  In  Figure  19. 
Each  call  has  a sensitivity  to  9 grey  levels  (this  may  be  changed  to 
achieve  desired  results),  each  level  described  In  Table  22. 

table  22.  Grey  Level  Seeponse  of  CCD 

Level  Voltage  Range  Symbol 

1 0.0  - 0.2 


2 0.2  - 0.3 

3 0.3  - 0.4  + 

4 0.4  - 0.5  * 

5 O.S  - 0.6  & 

6 0.6  - 0.7  X 

7 0.7  - 0.8  *» 

8 0.8  - 0.9  S 

9 0.9  - 999.9  ■ 


The  Image  la  clocked  from  the  array  with  any  desired  tranafer 
efficiency  or  contrast  ratld.  With  a transfer  efficiency  of  93X  end 
the  Interline-transfer  method  used,  the  observed  Image  becomes  that 
of  Figure  20.  Note  the  degradation  of  the  Image  as  well  as  the 
delay  encountered.  This  delay  Is  not  related  to  clock  frequency 
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other  than  through  the  transfer  efficiency.  When  pattern  tracking, 
to  maintain  no  delay  due  to  charge  diatributlon  changes,  the  trans- 
fer efficiency  oaist  be  well  above  99Z  as  illustrated  by  Figure  21. 

Of  course,  the  position  of  the  target  on  the  sensor  face  does  have 
sosw  bearing  on  the  target  pattern  degradation.  The  amount  of 
degradation  is  directly  proportional  to  the  number  of  transfers 
necessary  to  remove  the  image  scene  from  the  sensor.  The  simulation 
removes  the  image  by  clocking  the  cells  up  and  to  the  left.  There- 
fore, the  target  is  more  easily  found  in  the  upper  left-hand  comer 
than  the  lower  right-hand  comer. 

The  determination  of  the  minimum  transfer  efficiency  for 
target  recognition  begins  with  the  generation  of  a complicated 
background  on  which  to  view  the  target.  The  pseudo-random  back- 
ground of  Figure  22  provides  possibly  worst-case  results  for  pattern 
recognition.  The  output  of  the  program  using  the  target  shown  in 
Figure  23  is  tabulated  for  evaluation  in  Table  23.  The  contrast 
ratio  is  set  at  0 to  once  again  provide  wovst-case  resulti^i  and 
to  give  an  honest  feel  for  pattern  tracking  rather  than  positive 
contrast  tracking.  The  target  pattern  is  positioned  at  x ■ 50, 
y > SO,  resulting  in  approximately  100  transfers  for  each  cell. 

The  pattern  cell  tolerance  is  set  at  0.09  volts,  giving  approxi- 
mately one  grey  level  of  tolerance.  The  transfer  efficiency  must, 
therefore,  be  greater  than  99.90%  in  order  to  maintain  track  under 
these  conditions.  This  efficiency  is  of  course  «Icpendent  on  several 
factors  other  than  placsmsat  of  the  target  on  the  image  surface.  The 
complexity  of  the  target  Itself  is  an  Important  factor.  The  grey 
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level  tolerance  of  the  pattern  cell  can  be  raised  to  provide 
tracking  or  acquisition  at  lower  transfer  efficiencies.  However, 
the  possibility  of  locating  a "false"  target  Increases  with  an 
Increase  in  tolerance. 

An  estimation  of  the  lowest  transfer  efficiency  allowing 
acquisition  of  the  target  may  be  calculated  for  the  above  situation. 
Assuming  no  addition  from  preceding  charge  packets,  tolerance  of 
0.09  volts  and  pattern  location  of  x - 50,  y <■  50  ; the  largest 
number  of  transfers  is  106  and  the  value  of  that  cell  is  1 volt. 
Therefore, 

(n)“  (V)  - (V)  - (toi)  . (30) 


where  n is  the  transfer  efficiency,  N is  the  nu&ber  of  transfers, 
V is  the  value  stored  in  the  cell,  and  tol  is  the  cell  pattern 
tolerance.  This  gives 


(n)““  (i) 


.108 


1 - (0.09) 
0.91 


or. 


.108  iln  n 
iUi  n 

iln  n 


in  0.91 

ia  0.91 
108 

-0.00087 


n - 0.99913 


76 


.*•  ■*«  X*  MX**.X*  XX  mxt>*  *X.  MX  *S>4X.X>**;!*  IXX 

X «t  X4i  tx..x*»  •*  .•*  xt  X}  x4xxt  . •“ 

<4.uk  i t .xi:x«ii*x  t.xiixxx}xw  t^ktx'x.  •»  * !!!■  «.!??■ 

A I «x.x*xx*.ii  »x  »>aiNiM<.x<.aMxxx>'x  x<*»x  .«tXAX*»i.  • x**. 

t«>X«  «*«4k..«M  M «*tf<  •«,X4n«X*,  «.  iax.ll*  < ,*  *> 


+ 4X  XX*  X «4«*.** 
XX  ■ k X4X  X.XXtttXk  MW 

X xxx  *1  . ..tX'titxxii  4 d 

♦ ■4.x,  X*X..4tt  «.  _♦  .•X...X44fXX. 


X * A44XX.Xt*  .AX  WU  « ..»♦  MX.*  ‘X'.}  ••  X*li!*I*i*!'‘*SiL*2!y»li45li{<*S5*Ali!5  ixUl* 

«.  .lMUX4,x44k«iy4 1 ,x  ,4,y»  X.,ak<>Xy*444,«  k „4Ak»«*»Hyki»k  X 1 4 kt  xjijk 

Mk. .««xxki>4  4i4X.ki.  U «««k  *4Mx,i,.k  i<x.xk>xk4xxx^  *lil4l:!lx  Stx«x  x«x  xx*  xt, 

x.k.  XX  X4XX4.  Xl4. 
x'«k  4t».x.x*Sx..<'.il 
- • - • 4.4.X4 


. x.ts)  ji«X4k»<n}*  .Xkx  k 4 4a  i.«xH4X4  a.iti*  x*.x»4.4  jjtx.  bU.x  *,•  ■*  !JI*5  ! ’‘5!  StxxV 

.X44«X  X 4 XM44«i4.4kkXk.4  kk  4li  ,JS  4.  k 4444  X4  4iXk.  • ‘ f . J !!*  ?!  j!  v .jX44.X.»..XX.  kJA}, 

MAX*  X 4AkX4XX4«,X,X,A44k  .kAk*  . 4X  X4*  X..*i  *kX4l44  4 y , » 4 k 4 X 4 »«»■ JJX* X . . » J 
k *4  iy  ,X4kk  ..>.1  4 I ■kM4khk4tA  4,“»»  ■ 4 aM»  4 . AX  . « 4 XkX  < • tX  *44, XJ  . Xtk 

• 4,kykXM  ii  Xyy  4*4.  *4.11.  .iM4Xk  .*.XAM  4 , » . ik  44«,»X . y.  Xvkk  4i  k kk44  ••*.21'xi45'x7;i''”'4i«4"'’i 

y !(.  4 44  X 4k  *4  y xH  u i.ykJA  X4tk4  *4444  '>xX44.xk*tx»4l.4  •}SX  x»x  flxlx  is*!*  5 

x4..k.Xk^  .4x  . k . AkA,k.nkX4y»xx,4«x  . t«  ■ ! •iS********!l?**!*l!  ! yl***!!*4yJx»yi{Ixt  8 *"*■»  *♦ 
44Xi,Xjk4  4.<Ai4Hi,.  ^ kA  444  X »k*4A  4XA4.*4iX44  4 *t  Xk44  k 4 J J 4X  X 4 , , k X4  4XXXkk*  »■!  k 4*"* 
i4y.4tk  4 ,44yA4y*,iu  i4XkAt  .k  4kl4X,t  .4  k •(•ikl.  ,44k  xij  IJ’k’lJtkyyyik*!!  IIS*  . 

4 kkMX,  U 1 4*  X.  4X  44.  >»  «,  .44Vk4Xk4kMXXkX4klit4«  X X . t 4 « , 1 XX . XJ . 4 .4*  t k»4kk»k4X  XJ  JJX  ,, 
kiXX.  S4Xh  ■4X44„k.yA4A  Xm  »4  4 4 kX44A*k  A44444  4,  ♦X4XX  kXl'*  MJ  XkX  4.  JM  k4MX444X4kXXXXk4  ■ 
.X<UA  kAk.k.IlX  . kkX.44UI  4.  . kJX.l.,  & 44.  XXkl* * • X4 . . 4 X 4 4 4 . k4k  44X4Si»A  4X  ■■k44XX44»4 

4 k.  44. xi  .k  X •kXkv  4 4X4  4 X4X4X>.  *4444  k X 4 ».,»  44yX441  k4.5X444,JXX  ■•JXX  X 4,44  . kj  « 5 

• 4^44kXk4.X44XX  k4  444  X 4 ,4444.44  A.Xy.Xf  k4  44J  444  4.»4aS.;.4  4XXX.4.XklkXk.X..4  JMkkJ  .XIJXX* 

.4a  4X.t  ,X4H4kM*XXX.4  k 4i)y--.4..44«4Xk.44ky.. 44X4X44  *44j-‘X  •X44.X  4X4$k  4 tJX  XjXkkJ^  kXkJkXk 

k4.4kXX44  4X,k4  .4  y,.  ».(44  . k 4..4X»*4k44  k k 4X  XIX  *4  4kkk444  kk4(k4XX  . ■ J.i44k  4X4Xk44  4* 
a M X4Xk.l,4.  .44X44,  .,  4.4X44  4l4<<,  . X 4 XM  *4  444X4X  ««  I**  J w*!*!  i.  „ ly  I Pi^.^yllS  »!  ‘l  S 44 
X4ytk«k4|i.t44XX4.,4,  4*  4. 4X44444X44. 4XX4kX. 4a  .XXXX4X44X44X  XkX»  X ■4  4ll»XjX.4.4  XOMX*;  4 4._«_1_»4, 
k4«4.  X444«  4 X44/(  kk  4 4 4X44XXilA4  x#,  Xk.X4XX4k4  *X»»*,i;r.X  4ikJ  kj  *'f2XX«kk4X  •«X«XX4X.XXBX 

.X  4.  X44  4 44X  .k  XXX  kki.,4A4X4kXXk.kXtk44X4X44  .X4  kXX4|k  4X4  4 4JtXXXkX  .4k  Xl.k4.  jX  • 

.4ixx  X Xkyk.kMAX  .44  «,A44XX44XXk  kXMkXA  4 4 g . . XX  ••■..♦kJ'J  * • }i S! * * ? ! I , 5 * * t y *4 1 .}*  *y*4M41  ** 
k 44  X 4.4  k4  ..•uXA.  X A «4A<  4..k»k.44  ,ltxk.4  kik  kM.k  k»4lX4X.X  *S  xiA 

kk  Akk  ctkX44  a44A4Xk4XkX4k  4444  U 4 XX4  «4,.4X.S.4XkkX4Xkg4  4 41X444  »,*‘44*X  4 |l44.IIXk4  *4  t ■IBj 

.4-.t4  xk.Ak  44.4  ..kykl  4X  ,A>.44  ..4X4  X 4..4y  A. 44  4 XX44X4  k kk  Xk  ,^1X4.4  •„»4»k4IX4»»lx 

kkl*kk4k44X44k4X4  • . .X.  X,  44  4444AXk44X  <lt.  XX4  4Xk  4k4kk 

k 4X444k«*4»A4A4ak4  X4k  Ik  « 44  4444kX  kky..<.44XX  k44XAX4 
4 4 ,|Ak4k4  4yl...4  4kX444  X 4 xt.xx  Xt.  4/:4MX  4 AX4kX 

Ikl.  X 41  S44  k.xx  4XX4X4Ak  A .XX4I4X44 .kkt«  lk44kX4a44 

k ..SXk.XkX  « 4y.AA  .4  AlX  kNX4k4k44  X4.  ■..44k4.«kA  X 4,  ......  lyi  i t,  ii  i~y  iv 

ktXA.  k *t4*.X»k44,  tlk*i  Ak  4t,X4k  4MX4I 4X. . t 4l,k.  4XXka  4 {^ky*‘y*4ykIkXkXli 

a4.4kx,  I444XX444  n44XAXk«  k 44X  4X4k  .X.  .XkXXXt  XX  4 ‘X  4.kX4  , 44aXkAXi4Xk.4  4 4X  .k.4kJ4»X«XXX 
A y .xtx  Xklk  XI. X.  i,  4X4  4X1UI  444X4XA.X  ■ X X. ■..11X4414X4X44  I {Jiljfi 

4XtXX  44  kXA  Xk  44«4k4xXAtl  XkXSX.  .k.  . 4.4  .•XXk.X4X.44  4Xk  illik'i.il*  jJkSx*  !!“ 

4 Xm4  4 4.X.X  1.4,44  4 4X,kNkX  X.  4XX  XXXX.XX4.XX4. 4.  ■44XXX  k 4»  }4k  J,  k ;«4X»  *•♦*■19. 

. 44  . 4tx4l.Ai,Xit  ,44  X,  y44lli<  XX44.I4  XA»  X4  , 4 4.44X  X44  k_4X4Xt4k  ^ ,1 . y Si.  . . ?S*  ♦«  J! 

kX  kxX4X  XX4X  4 kA.XXXk.  At  4.  k«  Mi  4kl.X4t.j444  4X4M.  »»»  ♦ • i*  *?*i'.*  ll!**!!!lS*i**  i*i5y 

4.4,  4 X,k4  44444, 4A, 4 X44XX4  xk  X a 4X4  X ‘ ‘ - 

44a  4lk  iaiX4aa4yXA«4. 4k. 1114X44  . X.4tX4 

.••kaXkn  ,X44t.4<Xk  AX4ki.«  .X  4X  4X*  «44X. 
t A A kkkXX  k.kx  .i4X  Xx,ilXk4X4X4  14  k4Mt4 

4X4  ii44At44a4X  . kt4l4Xt,4ai  Nxk  k 44k4X.44»  yAyy-yjA,?  ♦ ^ .ji".' ' iayix  kx 

4A  Xa  alxx  1.44k  k 4kXX4t  kX  4 4 4,XtXX4  kXkatytX  XAX.4lk  . XIXX  JXjXX  X kkl  **4»X»  ♦ 4»»XXJ  If 

• X,..44«  ikakllM  .ttA4  XI  X..  4 ,kA4l4.  4..  4 X\X  XaX_,  4X  ..g  _Mk  4 4 X . . k 4 k 4k  4 k 44X»a  XHkI  •«  Xitkt 

-4.„  4.4  4 XAXylA.X  X 44...4-$X  X..,X44..  •-•"■‘Vli'J.il’kxIl*!  4i;l4«k  V.M 


A it  M •«  ♦M*  ♦ ♦ 

•♦aiikrSk.iiji  s 

. 1X4XX  X kkl  **4»X*  ♦ 4»IKXX.Sf 

r 4X  ..x'^xkWxVtk'kiki  » ‘J*;*  *?**!. li.*!*l* 

|X44iX  4k4S  X 4X  .lXt4  4X  t 

r kx  x.ki.kM.i  • 4k»ikik  .•.44a  xtix,xix<iat«| 

li* 


il‘:%ISl*tk!!A«tk5!jk  . kiix,.*..  a4.  4 •?«*::!•••>«  \K 


Figure  22.  CCD  Array  with  pseudo-random  background  and  target  Image. 
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TARGET  PATTERN  SIGNATURE 


i 


110  110  11 
10  1^^101 
0 1^11^10 
1^10  0 1^1 
1^10  0 1^1 
0 1 ^ 1 1 1 0 
10  1^^101 
110  110  11 

Figure  23.  8x8  CCD  Array  Used  aa  Target 
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Table  23.  Output  of  CCDSIM.SEA&CH 
Iteration  Output 

1 GOOD  SEABCU,  TRANSFER  EFFICIENCY  IS  100.00% 
THE  TARGET  WAS  LOCATED  AT  X « 50,  y - 50 

2 GOOD  SEARCH.  TRANSFER  EFFICIENCY  IS  99.90% 
THE  TARGET  HAS  LOCATED  AT  x - 50,  y » 50 

3 GOOD  SEARCH,  TRANSim  EFFICIENCY  IS  99.98% 
THE  TARGET  WAS  LOCATED  AT  x “ 50,  y - 50 

4 GOOD  SEARCH.  TRANSFER  EFFICIENCY  IS  99.97% 
THE  TARGET  WAS  LOCATED  AT  x - 50,  y - 50 

5 GOOD  SEARCH.  TRANSFER  EFFICIENCY  IS  99.96% 
THE  TARGET  WAS  LOCATED  AT  x - 50,  y - 50 

6 GOOD  SEARCH,  TRANSFER  EFFICIENCY  IS  99.95% 
THE  TARGET  WAS  LOCATED  AT  x - 50,  y - 50 

7 GOOD  SEARCH,  TRANSFER  EFFICIENCY  IS  99.94% 
THE  TARGET  WAS  LOCATED  AT  x » 50,  y « 50 

8 GOOD  SEARCH,  TRANSFER  EFFICIENCY  IS  99.93% 
THE  TARGET  WAS  LOCATED  AT  x • 50,  y - 50 

9 G00.D  SEARCH,  TRANSFETt  EFFICIENCY  IS  99.92% 
THE  TARGET  WAS  LOCATED  AT  x - 50,  y - 50 

10  GOOD  SEARCH,  TRANSFER  EFFICIENCY  tS  99.91% 
THE  TARGET  WAS  LOCATED  AT  x » 50,  y - 50 

11  RAD  CCD  SYSTEK—TARGET  LOST!  1! ! i I 
TRANSFER  EFFICIENCY  TOC  LOW  AT  99.90% 

LAST  AVAILABLE  POSITION  WAS  x * 50.  y - 50 
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This  corresponds  to  a transfer  efficiency  of  99,913X,  which  cor- 
responds very  closely  to  the  simulation  output. 

It  should  be  pointed  out  that  the  simulation  is  for  interline 
transfer.  If  line — or  Interline — transfer  Is  used,  the  number 
ent-ered  as  TREFF  is  actually  the  cube  of  the  transfer  efficiency 
for  a three-phase  device  since  three  transfers  arc  made  for  each 
sensing  element.  For  frame-  or  f leld- transfer , two  parts  must  be 
used.  For  the  transfer  to  the  storage  section  the  number  entered  is 
also  the  cube  of  the  transfer  efficiency  but  the  only  vertical 
transfer  Is  made,  so  half  as  many  transfers  are  made  for  a square 
device.  The  second  part  la  Identical  to  the  line  transfer,  having 
both  vertical  and  horizontal  transfer.  An  efficiency  of  at  least 
99.91Z  Is  presently  attainable  with  state-of-the-art  MOS  fabrication. 

For  actual  simulation  of  known  device  characteristics,  the 
contrast  of  target  to  background  may  be  varied  as  Is  necessary. 

Any  target  may  be  used  for  the  recognition  scheme.  The  only  pro- 
vision Is  that  the  target  approximation  muat  cover  a rectangular 
array  of  cells.  Digitized  target  patterns  are  easily  obtainable  and 
may  be  used  directly  by  the  scheme.  With  slight  revision,  the  pro- 
gram can  also  provide  target  acquisition  as  well  as  tracking. 

One  great  advantage  of  the  CCD  mechanlam  Is  that  it  is  self- 
scanning. This  mechanism  provides  the  ability  to  observe  the  cells 
one  at  a time  as  one  might  observe  bits  in  a shift  register.  The 
grey  levels  are  easily  chossn  and  implemanted.  The  uses  for  the 
CCD  Imager  as  a tracker  necessitate  compactness  and  agility.  In 
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order  to  process  the  Information  from  the  CCD  where  several  grey 
levels  are  not  necessary  the  information  may  be  encoded  directly  as 
it  Is  clocked  from  the  array  sensor.  If  the  Information  Is  bi- 
level (0  or  1),  as  shown  In  Figure  24,  the  black  and  white  informa- 
tion may  be  encoded  as  run-lengths  of  cells.  If  several  grey  levels 
are  necessary,  the  difference  between  grey  levels  is  encodable  with 
large  dynamic  sensitivity.  A great  deal  of  redundancy  exists  in 
black  and  white  imaging.  If  a priori  knowledge  is  available  on  the 
background,  a good  '’educated"  guess  can  be  made  as  to  the  probability 
of  various  run-lengths.  With  this  knowledge,  the  information  can  be 
encoded  with  a large  savings  in  the  number  of  bits  necessary  for 
reconstruction  for  the  image  scene.  A one-zero  representation  of 
Figure  24  is  presented  in  Figure  25.  If  this  "scene"  background  is 
assumed  to  have  a Poisson  probability  distribution  with  an  average 
run- length  of  5,  it  can  be  encoded  as  it  is  clocked  from  the  array 
using  the  Huffman  code  generated  from  the  appropriate  distribution. 
The  encoded  sequences  are  given  in  Table  24  for  lines  40  through  60 
of  Figure  25.  Since  the  target  is  not  accurately  represented  by  the 
Poisson  distribution,  it  is  easily  recognizable.  Even  though  the 
isuige  of  Figures  24  and  25  is  not  accurately  represented  by  the 
Poisson  distribution,  the  number  of  bits  necessary  for  description 
is  reduced  and  the  target  pattern  is  easily  found  in  the  sequence. 
Perhaps  a better  distribution  for  description  of  the  scene  is  the 
modified  binomial  distribution  previously  described.  A program  is 
presented  for  this  encoding,  along  with  the  code  to  be  used  in  its 
table  look-up  procedure  in  Appendix  B. 
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Figure  24.  CCD  Amy  vlth  bl-* level  representation  for  Huffman  Coding 
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Figure  25,  Oae-Zero  rep reeen ration  of  Figure  24  for  coding* 


Table  24.  Encoded  Sequences  for  Lines  40 
Through  60  of  Figure  25. 


1*011 

0*1010 

1*011 

0*0010 

0*1011000 

0*010 

1*011 

0*100 

1*1011000 

1*1010 

1*00110 

0*100 

1*000 

0*0010 

0*1010 

0*1011000 

1*1011000 

0*00111 

1*00111 

1*00111 

1*1010 

0*11 

1*00111 

0*1010 

0*00111 

0*1010 

1*1010 

0*011 

1*011 

1*100 

1*011 

0*011 

1*000 

0*00110 

0*011 

0*1010 

1*011 

0*010 

1*100 

1*00110 

1*010 

0*11 

1*00111 

0*100 

0*000 

0*010 

1*1011000 

0*1010 

1*100 

1*010 

1*1010 

0*00111 

1*00111 

0*011 

0*11 

0*0010 

1*011 

0*1011000 

1*0010 

1*1011000 

1*11 

0*11 

1*0010 

0*11 

0*00111 

0*1010 

1*1011000 

0*1010 

0*000 

1*011 

1*1011000 

0*100 

1*1010 

1*1011000 

0*000 

0*000 

1*00110 

0*010 

0*1011000 

1*000 

1*0010 

0*11 

1*11 

1*1011000  1 

( 0*1011000 

0*0010 

1*00111 

0*100 

0*011  1 

i 1*1011000 

1*1011000 

0*00110 

1*0010 

1*11  < 

0*00111 

0*0010 

1*0010 

0*1010 

0*00110 

1*1011000 

1*1010 

0*0010 

1*010 

1*1011000 

^0*1011000 

0*010 

1*100 

0*100 

0*011 

1*1011000 

1*00110 

0*00111 

1*1010 

1*1011000 

0*0010 

0*010 

1*00110 

0*0010 

C"’011 

1*010 

1*000 

0*00110 

1*1011000  i 

^ 1*00110 

0*000 

0*1011000 

1*0010 

0*010  \ 

0*1011000 

1*00111 

1*1010 

0*00110 

1*100  < 

[ 1*00111 

0*11 

0*0010 

1*000 

0*100 

1 0*1011000 

1*011 

1*1010 

0*1010 

1*11 

1^1*00111 

0*00111 

0*010 

1*00111 

0*100 

0*0010 

1*100 

1*1010 

0*1010 

1*0(J10 

1*100 

0*0010 

« 
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Table  24.  (Continued) 


1*000 


1*1011000 


1*010 
0*100 
1*000 
0*00110 
1*011 
0*1011000 
1*1011000 
0*1011000 
1*00111 
0*1011000 
. 1*1011000 


1*010 

0*00110 

1*010 

0*010 

1*1010 

0*00111 

1*000 

0*100 

1*00110 

0*000 

1*00110 

0*00111 


0*1011000 

1*1011000 

0*00111 

1*1011000 


0*1011000 

1*1011000 

0*1011000 

1*000 

0*000 


0*1010 

1*1010 

0*00111 

1*0010 

0*011 

1*1011000 

0*00110 


0*000 


0*011 

1*011 

0*1010 

1*011 

0*1011000 

1*1011000 

0*00111 

1*1011000 

0*1011000 

s. 

1*100 


1*010 

0*010 

1*011 


1*011 

0*00111 


1*000 

0*010 

1*00110 

0*100 

1*011 


1*000 
0*1011000 
1*1011000 
0*1011000 
1*00111 
0*1011000 
1*1011000 
0*00111 
1*1011000 
0*011 
1*1010 
0*00110 
1*1011000 
0*00110 
1*00110 
1*1010 
0*00110 
1*1010 
0*010 
1*100 
0*00111 
1*00111 
0*1011000 
1*1011000 
0*00  in 
1*1011000 


( 0*1011000 
I 1*00111 
0*000 
1*00110 
0*0010 


0*0010 

1*100 

0*1011000 

1*1011000 


1*0010 

1*010 

0*1011000 

1*010 

0*1011000 

1*1011000 

0*00110 

1*011 

0*0010 

1*00110 

0*00111 

1*00111 

0*10J1000 

1*00111 

0*1011000 

1*00111 

0*0010 

1*1010 

0*1011000 

1*00110 

0*1011000 

1*100 


0*00111 

1*100 

0*1011000 

1*0010 

0*010 

1*000 


1*000 

0*100 

1*011 

0*1010 

1*010 

0*00110 

1*011 


1*000 

0*0010 

1*1011000 

0*0010 

1*010 

0*100 

1*1011000 

0*000 


0*011 

1*1011000 

0*000 

1*1010 

0*010 

1*010 

0*011 


0*00110 
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Table  24  (Continued) 


1*0010 

0*1010 

1*010 

0*1011000 

1*010 

0*0010 

1*000 

0*000 

1*100 

0*00110 

1*1011000 

0*0010 

1*000 

0*00111 

1*11 

1*1011000 

1*100 

0*1010 

1*11 

0*0010 

1*1010 

0*1010 

1*11 

0*00110 

0*010 

1*00111 

0*00111 

1*00111 

0*11 

1*1011000 

0*00110 

1*000 

0*00111 


1*1011000 

0*11 

1*100 

0*000 

1*000 

0*100 

1*100 

0*100 

1*1011000 

0*00110 

1*00111 

0*00110 

1*00110 

0*100 

1*100 

0*011 

1*011 

0*1010 

1*1011000 

0*00111 

1*1010 

0*000 

1*00110 

0*1010 

1*100 

0*00111 

1*0010 

0*1011000 

1*00110 

0*0010 

1*00111 

0*0010 

1*1010 


0*100 

1*011 

0*011 

1*11 

0*011 

1*00110 

0*100 

1*0010 

0*1010 

1*1011000 

0*0010 

1*1011000 

0*011 

1*10110000 

0*1011000 

1*0010 

0*1011000 

1*00110 

0*00111 

1*1010 

0*1010 

1*00110 

0*1010 

1*00111 

0*11 

1*0010 

0*1011000 

1*100 

0*011 

1*11 

0*0010 

1*1011000 

0*010 


1*100 

0*0010 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  CCD  mechanlaa  lends  Itself  to  imaging  systems  that  are 
certainly  coaoensurate  with  the  device  capabilities.  Most  of  the 
parameters  of  the  SSID  as  used  as  sensors  in  trackers  are  not  the 
limiting  characteristics.  The  parameters  of  the  tracker  provide 
the  liodtations  in  most  cases.  System  degradation  due  to  sensor 
characteristics  are  mostly  charge  transfer  oriented.  This  study 
shows  the  CCD  to  be  capable  of  handling  up  to  10  MHz  operation 
which  is  sufficient  for  most  tracking  application. 

The  small  size  and  low  power  requirement  indicate  the  strong 
impact  on  image  tracking  where  limited  space  and  power  availability 
are  factors.  As  evidenced  by  this  study,  with  adequate  transfer 
efficiencies,  the  accuracy  of  the  position  of  the  target  to  which 
each  emergent  charge  packet  can  be  assigned,  sensitivity,  and  dynamic 
range,  give  these  sensors  wide  application  in  the  field  of  target 
acquisition  and  tracking. 

The  devices  are  easily  modeled,  giving  simulation  results  idiich 
can  be  used  to  predict  the  accuracy  of  a tracker  utilizing  the  CCD 
or  CID.  The  device  parameters  (transfer  efficiencies,  contrasts, 
target  patterns,  background  distribution),  ate  estimated  and  then 
used  by  the  simulation  program  to  predict  the  behavior  of  a tracker 
using  the  device  as  a sensor.  The  routine  is  easily  incorporated 
into  target  recognition  schemes  and  tracking  algorithms.  One  such 
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algorithm  is  the  correlation  method  for  target  or  image  tracking. 

The  study  of  this  tracking  technique  as  applied  to  the  CCD  simu- 
lation Is  suggested  as  a possible  extension  of  the  work. 

The  Issue  of  buried  channel  versus  surface  channel  mode  Is 
resolved.  Burled  channel  provides  several  major  advantages  in 
performance  and  also  almtllfles  device  design  and  operation.  The 
charge  transfer  efficiency  for  the  buried  channel  mode  Is  high  for 
the  full  range  of  electron  packet  size,  from  saturation  charge  of 
approximately  l0‘  electrons  to  the  order  of  10  electrons  or  less. 
Obviously,  this  provides  wide  dynamic  range  in  straightforward 
design  whereas,  with  the  surface  channel  mode.  It  is  necessary  to 
cope  with  the  level  of  the  variably  required  "fat  zero"  clumnel 
current.  One  criticism  of  the  burled  channel  device  is  that  the 
saturatlou  signal  charge  density  cannot  be  made  as  high  as  with  the 
surface  channel  device,  and  therefore  is  more  limiteo  in  dynamic 
range.  On  the  contrary,  the  buried  channel  device  offers  orders  of 
magnitude  higher  dynamic  range  by  virtue  of  the  relatively  low 
noise  levels. 

Image  blooming  can  only  be  remedied  by  the  Incorporation  of 
over-flow  sinks  which  must  be  applied  to  each  sensor  element  in 
order  to  accomplish  adequate  anti-blooming.  Even  when  an  image  is 
excessively  Intense,  smearing  occurs  in  device  types  where  charge  is 
transferred  through  illuminated  areas.  Both  frame- transfer  and 
line-transfer  device  organizations  have  this  problem;  the  interline- 
transfer  device  does  not  because  the  photoelemeut  sites  are  distinct 
from  the  transport  register.  Where  it  la  a problem,  the  device  must 


be  designed  and  operated  In  such  a way  that  the  smear-producing 
charge  transfer  process  is  carried  out  as  fast  as  It  Is  generated. 

A compromise  solution  to  the  problem  is  to  incorporate  only  one 
anti-blooming  charge  sink  per  sensor  column.  This  feature  can  pre- 
vent blooming  between  columns,  but  not  within  columns. 

The  recognition  scheme  used  in  this  work  is  a simple  pattern 
search.  With  a few  and Ificat ions,  the  recognition  scheme  could  be 
greatly  improved  to  encompass  more  sophisticated  pattern  searches 
utilizing  more  complex  digitized  target  patterns.  This  work  is 
meant  to  be  a foundation  for  the  simulation  of  the  device  itself  as 
well  as  its  incorporation  into  tracking  schemes. 

The  actual  calculation  of  the  parAmeters  read  into  CCDMAIN  can 
be  accomplished  to  further  the  simulation  before  the  sensor  is 
manufactured.  Determination  of  device  response  or  MTF  can  be  made 
from  design  specifications  before  the  prototype  is  developed. 

It  can  be  seen  that  for  a 100  x 100  CCD  sensor  array  which 
utilizes  10,000  photosensing  elements,  a 1 MUz  clockrate  will  yield 
a frame  rate  of  up  to  100  frame/sec.  This  freme  rate  Is  certainly 
high  enough  to  maintain  track  in  any  of  the  trackers  surveyed  in  the 
study.  The  transfer  efficiency  necessary  to  maintain  pattern 
recognition  in  the  simulated  scheme  (99.90%)  is  easily  attainable  at 
such  a clock  rate.  Therefore,  even  with  clock  rates  as  low  as  1 MHz, 
adequate  missile  crackers  utilizing  a CCD  sensor  are  realizable. 

When  utilizing  this  method  as  a missile  tracking  mechanism,  the 
need  for  data  reduction  is  evident.  The  scheme  for  encoding  the  run- 
lengths  of  white  and  black  cells  as  they  are  clocked  from  the  array 
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ia  Introduced.  The  Huffman  coding  scheme  is  used  with  the  table 
look-up  method  of  codeword  assignment.  Sample  encodings  are  made 
using  codes  derived  from  Poisson  probability  distributions  with  a 
mean  of  five  and  with  a modified  binomial  distribution  utilizing 
consecutive  runs  of  while  and  black  cells.  It  is  evident  that  more 
knowledge  of  the  background  scenery  encountered  is  desirable.  The 
more  closely  the  distribution  can  be  approximated,  the  more  reduction 
that  is  possible.  Other  distributions  should  be  investigated  and 
actual  digitized  scenes  should  be  analyzed  to  determine  the  optimum 
code  for  a particular  application. 

The  encoders  for  these  operations  are  very  easily  implemented 
using  simple  digital  electronics  and  could  be  Incorporated  into  the 
peripheral  electronics  of  the  clocking  waveform  generators  and  output 
amplifiers  of  the  device.  This  possibility  creates  a need  for 
special  chip  design  for  special  applications  such  as  trackers,  both 
with  onboard  tracking  and  onboard  transmission  of  data  to  ground 
stations. 

The  idea  is  already  being  used  for  storing  pages  of  literary 
information  by  encoding  the  outputs  of  linear  reticon  arrays,  these 
encoded  data  are  available  then  for  instant  recall.  Up  to  30  to  1 
reductions  of  such  information  are  attainable  in  present  systems. 

The  need  for  an  extension  of  this  work  into  the  area  described  above 
utilizing  charge-coupled  devices  and  area  arrays  is  apparent. 
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A Synopula  of  the  Coaputer  Kout inos 
Ueed  the  Paraaeter  Study 

All  of  the  software  developaent.  aiaulatlon,  data  recricval  aad 
data  reduc t Ion  are  done  on  the  Unlvac  1106  Multiprocessor  facility 
available  on  the  Mississippi  State  University  raapus.  A listing  of 
the  software  used  follows  this  synopsis.  All  of  the  prograas  except 
one  are  written  in  Fortran  IV  and  Fortran  V.  The  one  prograa  used 
for  data  retrieval  froa  the  9>track  IdM  coapatible  tape  is  written 
in  Unlvac  1106  assaabisr  language.  Since  the  data  tapes  use  16  bit 
«iorde,  the  word  lengths  are  changed  and  manipulated  to  achieve  the 
36  bit  word  length  necessary  for  use  on  the  Unlvac  facility. 

The  cooputer  software  are  discussed  in  three  separate  sections 
as  related  to  their  use  in  the  study.  The  data  retrieval  prograas 
and  subroutines  arc  discussed  first.  These  prograas  Include  file 
aanlpulation  on  tape,  word  length  translation,  error  calculation  and 
a Gould  plot  routine  for  displaying  either  the  data  or  the  error 
Incurred . 

The  next  section  listed  is  the  CCD  slaulatlon  using  a 100  x 100 
eleaent  array  and  facilities  for  elesMnt  aanlpulation  as  well  as  the 
variation  of  paraaeters. 

The  last  listing  Is  that  for  the  encoding  of  the  array  data  as 
It  is  clocked  froa  the  CCD.  The  encoder  uses  the  table  look-up 
procedure  and  includes  a listing  of  two  of  the  codes  used  for 
experlaental  evaluation  of  the  reduction  of  bit  sequence  length. 
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TAPE&EAD 

(Data  Retrieval  and  Evaluation) 

TAPEREAO. DRIVE! 

DRIVEI  la  uaed  to  apply  scale  or  calibration  factor  manipula- 
tion to  the  data  after  they  have  been  converted  to  the  correct  word 
length  and  stored  In  a buffer  (BUFF3).  The  errors  are  also  derived 
from  the  adjusted  tracker  output  and  target  position  for  each  time 
Interval.  Both  x-  and  y-position  errors  are  calculated  and  either 
may  be  transferred  to  the  plot  routine,  PLOTIT.  The  pregram  also 
contains  the  ability  to  test  for  tape  translation  errors. 

TAPEREAD.CONVKT 

Subroutine  BRKOUT  takes  the  36  bit  words  stored  In  the  960  word 
buffer  (BUFFI),  some  of  which  contain  more  than  one  niece  of  data 
information,  and  separates  them  out  Into  individual  computer  words 
using  the  Fortran  V FLD(I,J,H)  function.  The  Fi.n  function  allows 
bits  I through  J from  word  M to  be  transferr'id  to  a hit  sequence  of 
another  word.  These  %rards  are  stored  in  BUF.t3. 

TAPEREAD.ASMRD 

TREAD*  reads  the  data  tape  one  record  at  a time.  One  record  is 
equivalent  to  960  16  bit  words  of  426.67  36  bit  words.  Therefore, 
three  records  must  be  read  In  order  to  fill  an  Integer  number  of  36 
bit  words.  The  36  bit  words  are  stored  In  BUFFI  and  the  16  bit  words 


are  stored  In  BUFF2 
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COMVRT^  C«ke8  the  data  from  TREAD*  and  transforms  the  16  bit 
words  to  the  usable  36  bit  words  necessary  for  the  subroutine 
BRKOirr. 


TAPEREAD.PLOTIT 

Subroutine  PLOTIT  uses  the  data  transferred  from  DRIVEl  to 
generate  a Gould  plot  so  that  the  position  error  or  actual  position 
of  Che  target  can  be  evaluated  for  each  run  under  consideration. 

CCOSIN 

(CCD  Simulation) 

CC0S1M.CCI»«AIN 

CCDMAIN  Is  the  main  device  slauilation  program  on  which  any  Image 
can  either  be  generated  or  read.  The  parameters  for  the  device  are 
alec  read  into  the  routine^  affect  the  charge  transfer  and  target- 
backgound  contrast,  as  well  as  provide  the  grey  level  determination 
for  target  recognition.  Any  target  can  be  placed  on  the  Image  sur- 
face on  any  row-column  Iccatlon  desired. 

The  program  utilises  degradation  of  the  image  scene  as  deter- 
mined by  the  transfer  efficiency,  TREFF.  Both  charge  loss  and 
^•cldual  charge  addition  affect  the  clocked  array  image. 

CCDSIM.  SEARCH 

This  routine  searches  throughout  the  previously  clocked  array 
for  any  target  pattern  stored  in  the  target  array,  TARGET.  The 
location  of  the  target  area  is  given  and  the  transfer  efficiency  of 
the  device  at  which  the  array  was  clocked.  The  grey  levels  within 
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which  the  pattern  of  the  target  auat  lie  are  specified  to  enable 
the  u:;er  to  either  track  contraat  differences  or  specific  patterns. 
If  the  transfer  efficiency  is  sufTlciently  small  so  as  to  inhibit 
target  recognition,  the  CCD  is  said  to  have  lost  target  track.  In 
this  case,  the  laat  available  position  is  given  along  with  the 
transfer  efficiency  at  which  the  target  was  lost. 

CCDSIM.PATGEM 

PATGLM  generates  a random  sequence  of  run-lengths  of  white  and 
black  cells  used  in  the  evaluation  of  the  possibility  of  encoding 
the  array  data  aa  they  are  reauved.  Once  again,  the  erea  pattern  of 
the  target  is  placed  at  any  convenient  location  on  the  scene  image 
and  is  encoded  as  well. 


CCDSIM 

(Encoding) 

CCDSIM. ENCODE 

The  routine  used  for  encoding  the  scene  Image  utilizing  a sudi- 
fiad  binomial  probability  distribution  is  ENCODE.  The  run-length 
is  counted  and  identified  as  either  I's  or  O's  (black  cells  or  white 
cells).  The  neat  sequence  is  also  counted.  Therefore,  a run-length 
would  be  described  as  n I's  and  m O's  or  vice-versa.  The  probability 
is  assigned  aa 

P(S)  - p"q“  , (Bl) 


whete  p is  the  probability  of  a 1 occurring  and  q is  the  probability 
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of  • 0 occurrlag,  such  that  p q * 1.  In  the  run«  stade,  p * 0.6 
end  q • 0.4.  A run  scquMCc  might  appear  as 

11111000 

wd  characterized  aa  p^q^.  P(S)  is  thus  0.00498  and  is  associated 
vith  the  Huffman  codeword  00010110. 

CCDS1M.ENC0DE2 

If  the  run-lengths  are  to  be  characterized  by  either  a run  of 
I's  or  a run  of  O's^  but  not  both  consecutively,  a separate  encoder 
is  available.  ENC0DE2  counts  each  sequence  and  identifies  it  as  a 
run  of  either  I's  or  O's,  then  chooses  a Huffman  codeword  from  the 
code  library  to  transmit  or  store. 

CCDSIM.P0ISC0DE5 

POISCOOES  contains  the  codeword  library  used  in  the  lootc-up 
procedure  for  an  encoder  similar  to  EHC0DE2.  The  Huffman  codewords 
are  generated  from  a Poisson  probability  distribution  with  average 
wordlength  S.  This  codeword  library  is  taken  from  an  earlier  work 
in  the  area  of  redundancy  reduction. 

CCDSIM.HODBICODE 

This  library  contains  the  Huffman  codewords  generated  from  the 
modified  binomial  distribution  previously  described.  The  library 
will  be  given  in  App^dix  B.  This  codeword  listing  is  also  taken 
from  the  work  done  previously  in  redundancy  reduction. 
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CCOSIM 

(Pictorial  Description) 

CCOSIM. IMAGE 

The  next  prograa  and  subroutine  are  utilized  In  giving  a 
pictorial  description  of  the  CCD  Image  surface  before  or  after  the 
clocking  from  the  array.  The  picture  generation  shows  the  degra- 
dation of  the  IsMge  as  well  as  the  delay  produced  by  low  transfer 
efficiencies.  IMAGE  uses  the  data  on  the  CCD  to  generate  each  ''Ir*- 
ment  In  one  of  a possible  nine  grey  levels  that  may  be  specified  In 
any  Banner  so  desired. 


CCOSIM.  DATA! 

DAIAl  provides  the  array  description  to  be  used  by  IMAGE.  The 
data  are  read  f rosi  a temporary  storage  file  (having  previously  been 
cataloged  by  CCOSIM. CCDNAIN) . 


CCDSIM 

(Target  and  CCD  Characteristics) 

CCDSIM.  TARGET 

TARGET  is  the  element  used  to  describe  the  target  pattern 
utlllxed  54a  the  SEARCH  routine  as  well  as  In  the  main  p'^ogram  for 
storage  in  the  array. 


CCDSIM.CCDCHAR 


CCDCHAR  contains  the  initial  parameters  of  the  charge  coupled 
device.  These  are  used  for  characterizing  the  operation  of  the 
device  and,  In  addition,  target  placement,  contrast,  and  grey  level 
tolerances.  The  data  available  in  CCDCHAR  are  used  by  CCDSIM.CCDMAXN. 


TAPEBEAD  VARIABLE  LISTING 


BUFFI 

BUFF2 

BUFF3 

BUFF4 

BUFFS 

BUFF6 

BUFF? 

E&FZ 

ERFY 

SCLTM 

APPROX 

ERFXA 


BUFFER  STORAGE  FOR  THE  36  BIT  WORDS  CREATED  IN 
THE  SUBROUTINE  CONVRT*  FROM  THE  16  BIT  WORDS  IN 
BUPF2 

BUFFER  STORAGE  FOR  THE  16  BIT  WORDS  READ  FROM 
THE  DATA  TAPE 

DATA  INFORMATION  WORDS —CONTAINING  ALL  OF  THE 
INDIVIDUAL  TRACXER  DATA,  SUCH  AS  POSITION,  ZOTS 
DATA,  AND  TIME 

ZOTS  VOLTAGE  OUTPUT  INDICATION  OF  THE  X POSITION 
OF  THE  TARGET 

ZOTS  VOLTAGE  OUTPUT  INDICATION  OF  THE  Y POSITION 
OF  THE  TARGET 

TRACKER  INDICAIICW  OF  THE  TARGET  POSITION  IN  THE 
Y DIRECTION 

TRACXER  INDICATION  OF  THE  TARGET  POSITION  IN  THE 
X DIRECTION 

TRACKER-ZOTS  ERROR  IN  THE  X DIRECTION 
TRACKER-ZOTS  ERROR  IN  THE  Y DIRECTION 
TRACKING  TIME 

LINEAR  LEAST-SQUARES  FIT  TO  REMOVE  ELECTRONIC 
NOISE 

TRACKER-ZOTS  ERROR  WITH  ELECTRONIC  NOISE  REMOVED 


SCALE1,SCA1£2 

SCALE3,SCALE4 


SCALE  AND  CALIBRATION  FACTORS  FOR  ADJUSTING  BOTH 
TRACKER  POSITION  AND  ZOTS  INDICATION 


•* 
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TBEFF 

CCD 

XCCD 

TSTPT 

TKLESS 

CQNT& 

TOLER 

MAX 

NN 

XFLAG 

XSTPMl 

NXPOS 

HYPOS 

TSTPPl 


CCDSIM  VARIABLE  LISTING  USED  IN  TRANSFER 

TRANSFER  EFFICIENCY  ASSOCIATED  WITH  THE  CCD 
ASSUMING  ONLY  ONE  TRANSFER  PER  PHOTOSENSITIVE 
ELEMENT  (SEE  TEST  RESULTS  FOR  VARIATIONS  FOR 
DIFFERENT  TIMBERS  OF  PHASES) 

CCD  ARRAY  BEFORE  CLOOCING  THE  IMAGE  OUT 

CCD  ARRAY  INFORMATION  AFTER  CLOCKING  THE  IMAGE 
OUT 

DIGITIZED  TARGET  SIGNATURE 

AMOUNT  OF  REDUCTION  IN  TREFF  FOR  EACH  ITTERATION 
TARGET -BAOCGROUND  CONTRAST 

GREY’LEVEL  TOLERANCE  USED  IN  TARGET  RECOGNITION 

NUMBER  OF  ELEMENTS  SUBTENDED  BY  THE  TARGET 

NUMBER  OF  ELEMENTS  IN  THE  FIRST  ROW  OF  THE  TARGET 

TEST  FOR  RECOGNITION  OF  THE  TARGET 

LOWER  LIMIT  ON  THE  GREY -LEVEL  TOLERANCE  FOR 
TARGET  RECOGNITION 

X POSITION  OF  THE  TARGET 

Y POSITION  OF  THE  TARGET 

UPPER  LIMIT  ON  THE  GREY-LEVEL  TOLERANCE  FOR 
TARGET  RECOGNITION 


BESt  AVAIlABlt  COPY 


110 


IMlb  HUUlirtt  Illt  UtbKAUtU  IMAOL  t-UK 

A bivcri  lAHwtr  mITHAN  SHtClFlbU  CKbY  LbVLL  lOLEKAUCtS 

bObKOUUML  SLANCH(XCCU«1ELA0»MAX»NN»TSTP1  , I HCP'P-' » TOLCR ) 
UlMLAlbiOU  XbCiX  iUUt  lUUI  ftbIHTA  100)  tTSIKHK  loot*  1bTE>Hl(  lUO) 


0633luf=i7MAX 
TSTHMHlJT)  = Ti,»HT<  ljT)-IOLtH 
J:)  TSDH'iUJlisTblHUiuU^IOLbH 


20 

C»*« 

21 

c 

M- 

c 

L — 

1 

— « 

HAUt  lU  THb 


XT  unit— th-UlX-i;«  bQU(4)lt^  IllALlilt 
i«tXI  AULOXAtouT  lAHviCr\X£MtNI 


TKANtJ-^  U 


;»  iisiin 

JlSU 

2 J1SU1«1 

1<»0  10  2U 
J1=J1-J(42 


IfT3i  . iuu^4lNiiTo« 

iKtil.Ll.lliU)  VfW  TO  b 
20  •JTSOI 


00  TO  2 

La  ItiMtl 

lo(4^:AX=IA4•N 

lK(NN.AAA.t>T.MAX)  00  TO  0 
11=1 l«i 


UK  All  bcUMtllT  UOLb  UUT  COHKbbl^UI 4)  WiTTUN  OlVEN 
TOiXMAUCbb. 

UO  A Jl=Ul»10U 


>Jl=Jl-MOU(or  >1114)  A2 
ll=ll>>N«i 

a« 

« «1XH0;>=J1»M1/S-1 


!' 


besfavaiuble  copy 


prr.i 


uo  lt»J  iit=i»lbo  . 

lUl  WKikLllUl  (XCCUdHtUUI  •OK=1>10U1 

* i_  * 1 . . 


bO  KH'tUH*  •••yxu  ecu  btbftN'  - 1 AU6t  f LOST  1 1 1 1 ! ’ ) 

bl  i'Tlkku^Mt  t^HCifcNCt— -TOO  tOg  AT*  •2HF7.2* '** ) 

?OH.yAM*^*»»UAbV^HVAlUAyL£  HOblTlON  MAST'XS*  » !<*» bXf  • Ts*  • 14  J 
111  KttUKN 
trtu 


<*0uu  a(.AkUu  t>k'lClLt4CV  IS  iOU.UUX 

IHt  IAK(*Cl  WAt<  tOCAlLU  AT  X=  bU  Y=  bU 

.M)UU  SLAkLiif  IhAl.bftU  U-klLlUK?  IS-  99.VVi. 
IHL  iMkbc)  AkS  LUCAIlU  A1  X=  bU  Y=  bU 

«>OOU  U.AkC».»  YKAi.SKtk  LKFlcUriCY  IS  9^»96ti 
(Mfc.  lkk«H.I  *kS  LOCAIcU  AT  bU  Y=  bU 

.AtMM  bt-AiiUu  U-^-lClL^iCY  lb  W.97ii 

IHt  *mS  tuCAIcU  At  X=  bU  Y:  bU 

SA/OU  StAiiUti  1I.A<.bl-tk  tFFiCltNCV  IS  VS. 96b 
Ult  IAkbt.1  (.AS  LLCAItU  AT  A=  b(i  Yb  bU 

MtiA!  btAKCli*  IHAubFtK  U-riCltNtY  IS  VV.VbW 
IHt  IkKlM.)  WAS  tOCAIcLl  A I Ab  bO  V=  bU 

6V0U  StAkLil.  iHAnSI  Lk  tFl^  ICXtHCV  IS  VV.VHi 
•Nt  lAKVtt  kAS  tuCAItU  A1  bU  Y=  bU 

.UMAi  StAkLil.  IKA4.SI-tK  U’FiLiUH.Y  IS  VV.SSA 
IHt  lAHik-l  AkS  tOLAItU  AT  bU  Y=  bU 

MHAI  StAkLM.  lKA..S»tk  tl-klLitHCV  lb  W.V2A 
lilt  IkHOtl  mAS  LoLkIcU  Al  XS  bU  VS  bU 

-UMMJ  StAliLll.  ItAuiiALK  Ok  IL itlA. V IS  VV.VlS 
rilL  aAS  L^^LAIlU  Al  XS  bU  VS  bU 

UALi  CLU  SVbltH— fAkotf  i.Ubl!{!{! 
rNANSktk  OI  iLltl.LY>-IUU  tU»  Al  w.vas 
LASl  MVAltkmU.  I'usiuwi  »AS  kS  SO  vs  so 


t* 


CCD 

CODE 
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CCDSIM  VARIABLE  LISTING  USED  IN  ENCODING 


CCD  INFOSMATION  ARfiAY  TO  ENCODED 

CODEWORDS  ASSOCIATED  WITH  THE  VARIOUS  RUN-LENGTHS 
FOR  THE  J40DIFIED  BINOMIAL  DISTRIBUTION  OR  FOR 
THE  POISSON  DISTRIBUTION 

THE  MAXIMUM  NUMBER  OF  CONSECUTIVE  ONES  OR  ZEROS 
EXPECTED  IN  A RW 
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CTPTS 

NOPTS 

NGLEV 

NX 

NY 

TRAY 

DLxNE 


CCDSIM  VARIABLE  LISTING  USED  IN  IMAGES 

THE  OUTPOINTS  FOR  THE  DETERMINATION  OF  THE 
GREY  LEVEL  RANGE  OF  EACH  POINT 

NUMBER  OF  OVERPRINTS  USED  IN  THE  PRINTOUT 

NUMBER  OF  GREY  LEVELS  UTILIZED 

STARTING  X POSITION  OF  THE  ARRAY 

STARTING  Y POSITION  OF  THE  ARRAY 

THE  SYMBOLS  USED  FOR  EACH  OF  THE  OVERPRINTS 

DATA  POINTS  FOR  EACH  LINE  OF  THE  ARRAY 
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BCCumTV  clabmpication  or  tmib  rase 


WCUtMTV  Cl.MnriC»TeO««  or  TmU  Ot* 


The  purpose  of  this  study  Is  three-folil.  ‘i'he  first  part  consists  of 
evaluating  five  existing  trackers  which  utilize  solid-state  imagers  as 
sensors.  These  iaugers  are  of  three  types:  charge-coupled  devices  (CCDs), 
charge  Injiection  devices  (ClDs),  and  photodiode  arrays.  These  trackers  are 
investigated  as  the  contrast,  irradiance,  and  target  velocity  are  varied  on 
tlie  taage  plana. 

Since  the  charge-coupled  device  shows  great  potential  in  this  area,  the 
next  portion  of  the  study  is  devoted  to  the  evaluation  of  the  CCD  as  a 
sensor  in  a tracking  or  target  acquisition  algorithsi.  Parameters  are 
varied  in  a CCD  slssilation  schesM  to  provide  characteristics  for  the  device. 
Contrast,  target  pattern  signature,  transfer  efficiency,  and  grey  level 
element  tolerances  can  be  varied  for  the  sinulation.  The  simulation  pro- 
vides a 100  X 100  element  three-phase  interline  transfer  device. 

The  last  part  of  the  study  is  devoted  to  the  Investigation  of  coding 
the  video  data  as  it  is  clocked  from  the  array.  The  Huffun  coding  scheme 
is  used  to  eliminate  data  redundancy  by  run- length  encoding.  The  possi- 
bility of  target  acquisition  from  code  seqtaences  is  also  examined. 
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